4QA071498 


THEMIS 

T§  Of  18  058 


TEXAS  A&M  UNIVERSITY 

COLLEGE  STATION,  TEXAS  77843 


DISTRIBUTION  STATEMENT  A 
Approved  for  public  release; 
Distribution  Unlimited 


TEXAS  A&M  UNIVERSITY 
PROJECT  THEMIS 


Technical  Report  No.  55 


STATISTICAL  PERT:  AN  IMPROVED 
PROJECT  SCHEDULING  ALGORITHM 


Cynthia  S.  Dunn  and  Robert  L.  Sielken,  Jr. 


Texas  A&M  University 
Office  of  Naval  Research 
Contract  N00014-76-C-0038 
Project  NR047-700 


Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of 
the  United  States  Government. 

This  document  has  been  approved 
for  public  release  and  sale;  its 
distribution  is  unlimited. 


ATTACHMENT  I 


TecV  vn* C a \ Y v \ 


STATISTICAL ^EERT:  AN  IMPROVED 
S' PROJECT  SCHEDULING  ALGORITHM  « 


Cynthia  S./Dunn-aHt  Robert  L.  /Sielken,  Jr 


THEMIS  OPTIMIZATION  RESEARCH  PROGRAM 


INSTITUTE  OF  STATISTICS 
Texas  A&M  University 


Research  conducted  through  the 
Texas  A&M  Research  Foundation 
and  sponsored  by  the 

Q&t4ree^o&  Nyival  Research , 

//Contra<VNW^4-76-CH0J>38  / 

( /Lj^kct  HM47=7M — 


Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of 
the  United  States  Government. 

This  document  has  been  approved 
for  public  release  and  sale;  its 
dirtrlbutlon  is  unlimited. 


- ?X-  s s' 


ATTACHMENT  II 


3^7  3 VO 


A project  scheduling  algorithm  is  developed  and  illustrated.  For 


each  feasible  project  deadline  time  the  minimum  project  cost  and  corre 


sponding  optimal  deterministic  activity  durations  are  derived.  The  cost 


of  an  activity  is  assumed  to  be  a convex  piecewise  linear  function  of  its 


duration.  The  algorithm  is  based  upon  network-flow  techniques  including 


the  use  of  a labeling  procedure  which  preserves  complementary  slackness 


A computer  implementation  of  the  algorithm  is  documented 
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1.  INTRODUCTION 

This  paper  describes  a scheduling  algorithm  for  a project  composed 
of  "jobs"  or  "activities."  These  activities  are  represented  by  arcs  in 
a directed  network.  The  network  nodes  represent  events  in  time.  The 
activities  at  any  node  can  "commence"  as  soon  as  all  activities  "termin- 
ating" at  that  node  are  completed.  Associated  with  each  activity  is  an 
interval  of  possible  completion  times  and  an  associated  piecewise  linear 
cost  function.  Given  that  the  project  must  be  completed  by  a specified 
deadline  time,  the  algorithm  determines  the  individual  activity  completion 
times  which  minimize  the  total  project  cost.  Repeating  the  process  for 
all  feasible  deadline  times  yields  the  entire  project  cost  curve  and  asso- 
ciated optimal  activity  completion  times. 

For  example,  suppose  the  project  consists  of  activities  A,  B,  C,  D, 

E and  the  order  relations : 

A precedes  C and  D, 

B precedes  D, 

C and  D both  precede  E 

and  those  implied  by  transitivity.  The  corresponding  network  representa- 
tion is  shown  in  Figure  1 where  the  arcs  represent  activities  and  nodes 
are  events.  Notice  that  arc  F does  not  correspond  to  any  "real"  activity 
but  merely  represents  the  order  relation  that  A must  precede  D.  We  shall 
assume  that  such  dummy  activities  have  zero  completion  times  and  zero  costs. 

Using  this  network  representation  of  the  project,  the  problem  of  com- 
puting the  cost  curve  can  be  formulated  as  a network-flow  problem.  We  shall 
make  the  following  assumptions  about  the  network:  there  are  no  directed 
cycles,  and  each  arc  is  contained  in  some  directed  path  from  the  beginning 


± 
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node  (called  the  "source")  to  the  terminal  node  (called  the  "sink") . 

This  problem  can  also  be  formulated  as  a linear  programming  problem; 
however,  due  to  the  large  number  of  variables  and  constraints,  it  would  be 
impractical  storage-wise  to  solve  it  using  linear  programming  methods. 

D.  R.  Fulkerson  (1961)  has  formulated  a very  efficient  network-flow 
algorithm  for  solving  the  problem  with  a linear  activity  cost  function. 

In  this  paper,  Fulkerson's  algorithm  has  been  extended  to  accept  a convex 
piecewise  linear  cost  function  for  each  individual  activity. 


2.  PROBLEM  FORMULATION  AND  SOLUTION  PROCEDURE 


2.1.  Problem  Formulation 


The  cost  of  completing  an  activity  is  assumed  to  be  a convex  piecewise 
linear  function.  The  cost  curve  for  activity  I is  depicted  in  Figure  2. 
Note  that  the  allowable  completion  times  for  activity  I have  been  divided 

into  NK(I)-1  intervals:  [TIME(I,1),  TIME(I,2)],  [TIME(I,2),  TIME(I,3)] 

{TIME[I,NK(I)-2],  TIME [I ,NK(I)-1] } , {TIME[I,NK(I)-1] , TIME [I ,NK(I) ] } with 


TIME(I,  1)  <.  TIME(I,  2)  < ...  £ TIME  [I,  NK(I)]. 


(2.1) 


Here  we  interpret  TIME(I,1)  as  the  shortest  possible  completion  time  and 
TIME[I,NK(I) ] as  the  cheapest  completion  time.  Even  though  the  duration  of 
activity  I could  be  greater  than  TIME[I,NK(I)] , such  durations  would  be 
needlessly  expensive  and  hence  TIME[I,NK(I)]  is  the  practical  upper  bound 
on  the  duration  of  activity  I.  The  intermediate  times,  TIME(I,2),  ..., 

TIME [I ,NK(I)-1] , will  be  called  "breakpoints." 
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Breakpoints  arise  when  there  are  alternative  methods  of  performing  an 
activity.  These  methods  do  not  differ  In  the  end  result,  but  they  do  differ 
in  the  amount  of  time  they  take  and  thei^  cost.  For  example,  suppose  that 
snow  plows  rent  for  a fixed  $200/day  and  cost  a varying  amount  per  hour  to 
operate  depending  upon  the  speed  at  which  they  are  operated.  A correspond- 
ing activity  cost  curve  might  be  as  in  Figure  3 where  the  "breakpoints" 
correspond  to  the  use  of  different  numbers  of  plows. 

The  cost  for  completing  activity  I in  time  TIME(1,M)  is  COST(I,M) 
which  satisfies 

C0ST(I,1)  > COST (I, 2)  > ...  > COST[I,NK(I) ] . (2.2) 

Furthermore,  letting  C(I,H)  represent  the  rate  of  decrease  in  the  cost  of 
activity  I on  the  M**1  interval  implies 

cn  n - cost(i.i)  - cost (i ,2) 
v * ' TIME  (1, 2)  - TIME (1 ,1)  ’ (2.3) 


C(I,NK(I)-1] 


COSTf I.NK(I)-1]  - COSTf I ,NK(I) 1 
TIME[I,NK(I) ] - TIME[I,NK(I)-1]  * 


The  convexity  of  the  piecewise  linear  cost  function  implies  that 


C(I,NK(I)-1]  < C[I,NK(I)-2]  < ...  £C(I,1).  (2.3a) 

Let  XACT(I)  represent  the  duration  time  for  activity  I.  This  duration 
time,  XACT(I),  can  be  decomposed  as 


NK(I)-1 

XACT(I)  - l XACT(I,11) 

M-l 


(2.4) 


where 


XACT(I,1)  - min [TIME (I, 2) , XACT (I) ] 


XACT(I,M)  = min{TIME(I,M+l)  - TIME (I, M) 


max[0,XACT(I)  - TIME(I,M)]} 


For  example,  suppose  that  in  Figure  4 XACT (I)  =25,  then 


XACT (1,1)  = min[TIME(I,2) , XACT(I)] 


min [ 10 , 25] 


XACT (1,2)  = min{TIME(I,3)  - TIME(I,  2),  max[0,  XACT(I)  - TIME(I,2)]} 


min{20  - 10,  max[0,  25  - 10]} 


XACT (I, 3)  = min{TIME(I,4)  - TIME(I,  3),  max[0,  XACT(l)  - TIME(I,3)]} 


min{30  - 20,  max[0,  25  - 20]} 


XACT (I)  = l XACT(I,M) 
M=1 
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The  total  cost  associated  with  duration  time  XACT(l)  for  activity  I 


NK(I)-1 

KK(I)  - l C (I , M) XACT(I ,M)  (2.7) 

M=1 

where 

KK(I)  - COST (I ,1)  + C (1,1) TIME (1,1).  (2.8) 


The  total  project  cost  is 

NK(I)-1 

l [KK(I)  - l C ( I ,M) XACT ( I , M) ] . (2.9) 

I M=1 


Let  the  node  time  XNODE(K)  be  the  "length"  of  the  longest  path  from 
the  source  node  to  node  K when  the  "length"  of  an  arc  (activity)  is  its 
completion  time.  Thus,  for  example,  in  Figure  1 


XNODE(l)  = 0, 

XNODE(2)  = A, 

XN0DE(3)  = max(B,  A + F) , 

XN0DE(4)  * max (A  + C,  B + D,  A + F + D) , and 
XN0DE(5)  = max (A  +C+E,  B+D+E,  A+F+D+E). 

If  activity  I originates  at  node  Oj,  terminates  at  node  T , and  takes 
XACT(I)  units  of  time,  feasibility  requires  that 

XNODE(Oj)  + XACT  (I)  < XNODE^).  (2.10) 

Note  that  the  time  to  complete  the  entire  project  is 

XNODE(SINK)  - XN0DE( SOURCE) . 


I I 


u 


In  what  follows 


XNODE( SOURCE)  = 0 


without  loss  of  generality 


The  problem  is  to  minimize  the  total  project  cost  (2.9)  subject  to  the 


condition  that  the  project  is  completed  by  a specified  time  LAMBDA.  This 


C(I,M)XACT(I,M)]} 


min{PCOST(LAMBDA) 


subject  to  the  constraints 

NK(I)-1 

XNODE(Ot)  + I 
1 M=1 


XACT ( I , M)  - XNODE(T  ) < 0,  all  I 


XNODE(SINK)  < LAMBDA 


XACT (I, M)  £ 


XACT (I, M)  _>  L(I,M) 


TIME(I,2)  M - 1 

TIME(I,M+1)  - TIME(I,H)  M = 2 


TIME (1,1) 


0 - the  origin  node  of  activity  I 


T » the  terminal  node  of  activity  I 


Since  the  addition  or  subtraction  of  a constant  ih  the  objective  function 


does  not  change  the  problem,  we  can  represent  the  objective  function  as 


NK(I)-1 


max£  l C ( I , M)  XACT ( T , M) . (2.19) 

I M=1 


We  shall  solve  this  problem  for  all  feasible  values  of  LAMBDA.  The 
minimum  feasible  value  of  LAMBDA,  LMIN,  is  the  length  of  the  longest  path 
from  the  source  to  the  sink  when  the  XACT(I)'s  are  at  their  lower  bounds, 
XACT(I)  = TIME(I , 1)  for  all  I.  The  maximum  value  of  interest  for  LAMBDA, 
LMAX,  is  the  length  of  the  longest  path  from  the  source  to  the  sink  when 
the  XACT ( I ) ' s represent  the  cheapest  practical  times,  XACT(I)  « TIME[I,NK(I) ] 
for  all  I.  Thus,  for  a given  LAMBDA  such  that 


LMIN  < LAMBDA  _<  LMAX, 


the  constraints  (2.13)  - (2.16)  are  feasible.  The  proof  for  this  and  all 
other  underlying  theorems  presented  in  the  problem  formulation  and  algorithm 
are  found  in  Chapter  3,  Section  2.  We  shall  refer  to  the  problem  given  in 
(2.13)  - (2.19)  as  the  Primal  Problem. 

In  the  Primal  Problem,  dummy  activities  may  be  assumed  to  have  times 
and  costs  equal  to  zero. 

’ 2.2.  The  Dual  Problem 

The  standard  duality  theory  for  linear  programming  implies  that,  if 
the  primal  problem  has  the  form 

T 

max  c x 

subject  to  the  constraints 


Ax  £ b , 


(2.20) 


then  the  corresponding  dual  problem  is 


T 

min  b w 

subject  to  the  constraints 
*T 

A w * c 

w>0,  (2.21) 

see  for  example  Hadley  (1962) . Writing  our  Primal  Problem  in  the  form 
(2.20)  implies  that  our  dual  problem  can  be  written  as 


mini  LAMBDA  • V + ][]»U(I,M) 
IM 


G(I,M)  - XlL(I,H)  . H(I,M)  ] 
IM 


(2.22) 


subject  to  the  constraints 


F(I)  + G(I,M)  - H(I,M)  * C(I,M) 


all  I,M 


(2.23) 


l F(I)  - l F(I) 


IeO-^K 


IeTT=K 


K = node  j SOURCE, SINK 
K = SINK, 

(2.24) 


F(I),  V,  G(I,M),  H(I,M)  > 0. 


(2.25) 


Note  that  the  coefficients  in  (2.21)  of  the  s^  dual  variable  are  the  co- 
efficients in  the  sC^  primal  constraint,  so  that,  there  is  a natural  one- 
to-one  correspondence  between  primal  constraints  and  dual  variables. 

The  dual  problem  (2.22)  - (2.25)  can  be  interpreted  as  a flow  problem 
for  the  project  network.  The  dual  variable,  F(I) , associated  with  constraint 
(2.13)  is  the  flow  for  the  I activity.  The  constraint  (2.24)  implies 
that  except  for  the  source  and  the  sink  the  flow  going  into  a node  equals 
the  flow  coming  out  of  that  node.  Thus  at  all  nodes  other  than  the  source 
and  the  sink  we  have  conservation  of  flow.  The  total  flow  of  the  network  is 


SINK 


and  V is  the  dual  variable  associated  with  constraint  (2.14)  for  a fixed 


LAMBDA 


The  dual  variables  G(I,M)  and  H(I,M)  are  associated  with  the  upper 


and  lower  bounds  for  XACT(I,M)  respectively 


Rearranging  (2.23),  we  have  an  equation  of  the  form 


For  a fixed  value  of  f,  we  have  c - f = r 


In  (2.22)  we  want  to  minimize  an  expression  of  the  form 


or  equivalently  using  (2.27) 


Since  g = h + r and  both  g _>  0 and  h ^ 0,  making  h as  small  as  possible 


implies 


Correspondingly 


Thus 


g ■ max(0,  c - f) 


h * max(0,  f - c) 
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and  correspondingly 


G(I,M)  =■  max[0,  C(I,H)  - F(I)  ] , 


(2.28) 


H(I,M)  = max[0,  F(I)  - C(I,M)]. 


(2.29) 


Using  (2.28)  and  (2.29)  the  dual  becomes 


win (LAMBDA  • V + ££u(I,M)  • max[0,  C(I,M)  - F(I)] 
IM 


- JjLtt.10  • max[0,  F(I) 
IM 


- C(I , M)  ] } 


(2.30) 


subject  to  the  constraints 


l F(I)  - l F(I) 


IeOj'K 


ieTj-K 


■{: 


K = node  y SOURCE,  SINK 


K - SINK, 


F(I) , V > 0. 


A key  observation  at  this  point  is  that  for  all  (I,M) 


U(I,M)max[0,  C(I,M)  - F(I) ] - L(I,M)max[0,  F(I)  - C(I,M)] 

(2.31) 


is  a convex  piecewise  linear  function  of  F(I)  as  sketched  in  Figure  5. 


The  convexity  of  (2.31)  follows  from  U(I,M)  ^L(I,M).  Furthermore,  since 


the  sum  of  convex  piecewise  linear  functions  is  also  a convex  piecewise 


linear  function,  it  follows  that 


£U(I,M)max[0,  C(I,M)  - F(I) ] - ^L(I,M)max[0,  F(I)  - C(I,K)] 

M H (2.32) 


IE 
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is  a convex  piecewise  linear  function  of  F(I)  as  sketched  in  Figure  6. 

The  piecewise  linear  behavior  of  (2.32)  suggests  the  following  de- 
composition of  F(I): 

NK(I) 

F(I)  - I F(1,J)  (2.33) 

J-l 

where 

0 < F(I,J)  < C[I,  NK(I)  - J]  - C(I,  NK(I)  - J + 1]  (2.34) 

and 

C(I,0)  = » 

C[I,  NK(I)  ] =0. 

(also  see  Figure  7). 

Using  (2.33)  and  the  definition  of  L(I,M)  in  (2.18)  the  dual  objective 
function  becomes 


IIU(I.M)  • C(I , M) 

IM 

P NK(I)-M  1 

£U(I,M)  I l F(I , J) 

- TIME (I , 1)  • F[I,  NK(I) ) 

M L J-l  J 

i 

Using  the  definition  of  U(I,M)  in  (2.17),  the  dual  objective  function  can 
be  further  simplified  to 

LAMBDA  • V + ^U(I.M)  • C(I,M)  - 2£TIME[l,NK(I)+l-j]F(I,J) . (2.35) 

1M  1J 

Since  JJU(I,M)C(I,M)  is  a constant,  (2.35)  is  equivalent  to 
IM 


l 
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LAMBDA  • V - £JTIME(I,NK(I)  + 1 - J)F(I,J). 
IJ 

The  final  form  of  the  dual  problem  Is 


(2.36) 


min (LAMBDA  • V - ££TIME(I ,NK(I)  + 1 - J)F(I,J)} 
IJ 


subject  to  the  constraints 


i rjFd.j)!-  i r jp(i.j)"i  =| 

}i*K  u J ieTi=K  L J J t-v  k 


(2.37) 


- node^ SOURCE, SINK 


(2.38) 


0 _<  F(I, J)  je  C[I,  NK(I)  - J]  - C(I,  NK(I)  - J + 1]  (2.39) 


0 < V. 


(2.40) 


This  dual  problem  can  be  solved  as  a network  flow  problem.  The  origi- 
nal project  network  is  enlarged  by  adding  one  arc,  say  ARC(I,J),  for  each 
F(I,J)  so  that  the  project  network  has  NK(I)  arcs  from  0^.  to  T correspond- 
ing to  each  activity  I.  The  flow  F(I,J)  through  the  ARC(I,J)  has  the  cap- 
acity restriction  indicated  in  (2.39).  The  dual  problem  is  now  solved  by 
constructing  a network  flow  from  the  source  to  the  sink  that  minimizes 
(2.36)  subject  to  the  capacity  restriction  in  (2.39).  Refer  to  Figure  8 
for  an  example. 


2.3.  The  Network-Flow  Algorithm 


An  efficient  algorithm  for  generating  the  entire  project  cost  curve 
is  given  below.  The  algorithm  is  initially  sketched  in  general  terms  and 
then  presented  in  detail. 


1 


2.3.1.  The  Sketch 


The  algorithm  starts  with  the  largest  LAMBDA  of  interest,  LMAX,  and 


sequentially  determines  the  LAMBDA  corresponding  to  each  breakpoint  of  the 


convex  piecewise  linear  project  cost  function  PCOST (LAMBDA)  defined  in 


(2.12).  Corresponding  to  each  of  these  LAMBDA  values,  node  numbers  XNODE(K) 


are  produced  representing  optimal  event  times  in  the  project  network.  The 


corresponding  optimal  activity  times  XACT(I)  for  the  project  are 


XACT(I)  = min{TIME [I,  NK(I) ] , XNODE(T  ) - XNODE(0T)>  (2.41) 


since  the  activity  cost  is  a non-increasing  function  of  the  activity  time 


Recall  that  the  proofs  of  claims  are  given  altogether  in  Chapter  3 


The  algorithm  is  based  on  a labeling  process  in  which  labels  are  as 


signed  to  some  of  the  nodes.  In  general,  the  procedure  is  a systematic 


search  for  a path  from  the  source  to  the  sink  having  certain  desired  pro' 


per ties.  Flow  along  this  path  may  travel  through  arcs  either  in  the  same 


direction  as  their  orientation  or  in  the  opposite  direction.  Such  flows 


will  be  called  forward  and  reverse  flows  respectively.  Roughly  speaking 


a reverse  flow  is  really  only  a reversal  or  re-routing  of  earlier  flow  in 


the  forward  direction.  No  net  flow  in  the  reverse  direction  is  allowed 


The  labeling  process  is  started  with  a feasible  and  optimal  solution 


to  the  primal  and  dual  problems  for  LAMBDA  = LMAX.  The  initial  node  times 


are  found  by  setting  the  activity  times  equal  to  their  upper  bounds.  These 


initial  XNODE(K)'s  and  the  initial  flow  - F(I,J)  = 0 for  all  (I.J) 


satisfy  the  following  properties 


J- 


I I 

g 


e 
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ABAR(I,J)  < 0 ’•> F(I  ,J)  « 0,  and 


(2.42) 


ABAR(I,J)  > 0 ■>  F(I , J)  = C[I,  NK(T ) - .1 1 - C[I,  NK  ( I ) - .1  + 1 

(2.43) 


where 


' : 


ABAR  ( I , J)  = TIME  [I,  NK(I)  + 1 - J]  + XNODE(0I)  - XNODE(Tj)  . 

(2.44) 


Note  that  no  restrictions  are  placed  on  F(I,J)  when  ABAR(I ,J)  “ 0.  Hence- 
forth, the  properties  (2.42)  and  (2.43)  will  be  referred  to  as  the  "optima- 
lity properties"  for  LAMBDA  = XNODE(SINK) . These  optimality  properties 
imply  that  complementary  slackness  holds  and  that  the  flow  F(I,J)  minimizes 
(2.36). 

The  labeling  process  has  been  divided  into  two  parts  called  the  first 
and  second  labelings,  respectively.  In  both  of  these  procedures  we  have 
freedom  to  label  with  respect  to  complementary  slackness  since  we  work 
exclusively  with  arcs  having  ABAR(I,J)  = 0.  The  first  labeling  seeks  a 
path  from  the  source  to  the  sink  composed  of  infinite  capacity  arcs,  i.e. 
those  corresponding  to  J = NK(I).  If  such  a path  is  found,  the  algorithm 
terminates  since  the  Primal  Problem  will  be  infeasible  if  the  current  value 
of  LAMBDA  is  decreased.  If  no  such  path  is  found,  we  go  on  to  the  second 
labeling  in  which  we  search  for  a path  from  the  source  to  the  sink  having 
the  following  desired  properties:  for  all  forward  arcs  of  the  path  ABAR(I,J) 
= 0 and  F(I,J)  is  less  than  its  upper  bound  in  (2.39);  for  all  reverse  arcs 


1 


II 


of  the  path  ABAR(I,J)  = 0 and  F(I,J)  > 0.  If  at  the  end  of  the  second  label- 
ing the  sink  has  been  labeled,  we  say  "breakthrough"  has  occurred. 
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If  breakthrough  occurs,  then  the  minimum  arc  capacity  along  the  path  is 
determined,  say  CAP (S INK) . The  old  flow  F(I,J)  is  changed  by  adding 
CAP (SINK)  to  the  amount  of  all  forward  flows  on  the  path  and  by  subtract- 
ing CAP (SINK)  from  the  amount  of  all  reverse  flows  on  the  path.  This 
new  flow  still  satisfies  the  optimality  properties  and  is  interpreted  as 
an  alternate  optimal  dual  solution  for  the  current  LAMBDA® XNODE (SINK) . 

On  the  other  hand,  if  the  sink  has  not  been  labeled  at  the  end  of  the 
second  labeling,  we  say  "nonbreakthrough"  has  occurred.  When  this  happens, 
the  old  dual  variables  are  optimal  for  the  old  primal  problem  and  no  new 
alternate  dual  solution  can  be  found.  In  this  case  the  node  numbers 
XNODE(K)'s  are  changed  by  subtracting  a positive  quantity  DEL  from  all 
XNODE(K)  corresponding  to  unlabeled  nodes  K.  This  does  not  change 
XNODE (SOURCE)  = 0 but  reduces  XNODE(SINK)  = LAMBDA  by  DEL.  Through  (2.41), 
these  new  node  times  imply  a set  of  optimal  activity  times  for  the  new 
LAMBDA  where 

new  LAMBDA  = old  LAMBDA  - DEL. 

The  definition  of  DEL  guarantees  that  the  new  XNODE(K)'s  and  the  old 
F(I,J)’s  still  satisfy  the  optimality  properties.  Hence,  when  nonbreak- 
through occurs,  we  have  identified  the  point  on  the  project  cost  curve 
corresponding  to  the  new  LAMBDA. 

The  second  labeling  can  terminate  only  in  breakthrough  or  nonbreak- 
through. After  either  of  these,  the  entire  labeling  process  is  repeated. 

. *j  j 

2.3.2.  The  Details 

Initially,  the  algorithm  sets  each  activity  time  to  its  smallest 
(most  expensive)  feasible  value  and  determines  the  corresponding 


minimum  feasible  project  completion  time  (deadline  time  l,MTN) . Then 


the  algorithm  sets  each  activity  completion  time  to  its  largest  (cheapest) 
feasible  value  and  determines  the  corresponding  minimum  project  cost  and 
maximum  completion  time  of  interest  (deadline  time  LMAX) . 


The  iterative  procedure  is  begun  with  the  node  times  XNODE(K)  cor- 
responding to  all  activity  completion  times  at  their  largest  (cheapest) 
values  and  all  flows  F(I,J)  equal  to  zero.  These  node  times  and  tlows 


satisfy  the  optimality  properties 


A.  Labeling  Process.  During  this  routine,  a node  is  considered 


to  be  in  one  of  three  states:  unlabeled,  labeled  and  unscanned,  or 


labeled  and  scanned.  Initially  all  nodes  are  unlabeled 


In  general,  a node  label  has  four  parts  [A,  B,  C,  D]  when  the  node  is 


being  labeled  because  it  is  at  "the  other  end"  of  an  arc  associated  with 


If  "the  other  end"  is  the  terminal  node  T_,  then  the  label 


contains 


A=0_,  B=J,  D*  maximum  allowable  flow,  and 


C m 0 [denoting  that  flow  will  be  in  the  forward  direction 


If  "the  other  end"  is  the  origin  0_,  then  the  label  contains  A * T 


J,  D **  maximum  allowable  flow,  and  C ■ 1 [denoting  that  flow  will 


be  in  the  reverse  direction  (T 


1.  First  Labelin; 


Assign  the  source  node  the  label  [ 


CAP (SOURCE) 


In  general,  select  any  labeled,  unscanned  node,  say 


node  n,  and  search  for  all  unlabeled  nodes  Tt  such  that  n * 0_  and 


ARC[I,  NK(I) ] is  an  arc  with 


ABAR[I,  NK ( I ) ] = 0. 


(2.45) 


Label  such  nodes  with  [0^,  NK(I),  0,  CAP(Tj)  * «].  Such  Tj's  are  now 
labeled  and  unscanned,  and  node  n Is  labeled  and  scanned.  Repeat  this 
step  until  either  the  sink  node  Is  labeled  and  unscanned,  or  no  more 
nodes  can  be  labeled  and  the  sink  node  Is  unlabeled.  In  the  former  case, 
terminate  the  algorithm.  In  the  latter  case,  go  on  to  the  Second  Label- 
ing. 

2.  Second  Labeling.  Nodes  that  were  labeled  from  the  First  Label- 
ing retain  their  labels.  However,  all  nodes  revert  back  to  an  unscanned 
state.  The  general  step  is  to  select  any  labeled,  unscanned  node,  say  n. 

(1)  Scan  n for  all  unlabeled  nodes  Tj  such  that  n * 0j.  For 
each  such  node  find  the  J (if  one  exists)  such  that  both 

ABAR(I,J)  = 0 (2.46) 

and 

F(I, J)  < C[I,  NK(I)  - J]  - C[I,  NK(I)  - J + 1];  (2.47) 

then  assign  node  the  label  [0j,  J,  0,  CAP(Tj)]  where 

CAP(T  ) =»  mint  CAP  (Oj) , C[I,  NK(I)  - J]  - C[I,  NK(I)  - J + 1]  - F(I,J)} 

(2.48) 

so  that  is  now  labeled  and  unscanned.  If  no  such  J exists,  the  node 
Tj  is  not  labeled. 

(ii)  Scan  n for  all  unlabeled  nodes  0^  such  that  n * T^.  For 
each  such  node  0^  find  the  J (ii  one  exists)  such  that  both 


ABAR(J , J)  = 0 


(2.49) 
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and 

F(I,J)  > 0;  (2.50) 

then  assign  node  0j  the  label  [Tj,  J,  1,  CAP^)]  where 

CAPtOj)  - mlntCAPC^),  F(I,J)]  (2.51) 

so  that  0^  is  now  labeled  and  unscanned.  If  no  such  J exists,  the  node 
0T  is  not  labeled. 

Repeat  the  general  step  until  either  the  sink  node  is  labeled  and 
unscanned  (breakthrough) , or  no  more  nodes  can  be  labeled  and  the  sink 
node  is  unlabeled  (nonbreakthrough).  If  breakthrough  occurs,  go  on  to 
routine  B;  if  nonbreakthrough  occurs,  go  to  routine  C. 

B.  Flow  Change.  The  labeling  process  has  resulted  in  breakthrough. 
The  sink  node  will  have  a label  of  the  form  [0  , J,  0,  CAP (SINK)].  The 
total  network  flow  will  now  be  increased  by  CAP(SINK).  The  flows  are  up- 
dated as  follows.  Add  CAP(SINK)  to  F(I,J);  then  go  on  to  node  n = 0^ 
and  its  label.  The  general  step  for  node  n depends  on  its  label  and  is: 

1.  Label  - [Oj,  J,  0,  CAPO^)].  Add  CAP(SINK)  to  F(I,J) 
since  this  additional  flow  along  ARC(I,J)  will  be  for- 
ward flow  from  0^  to  n * T^.  The  next  node  to  consider 
is  n ■ 0^. 

2.  Label  - [Tj,  J,  1,  CAP(0T)].  Subtract  CAP (SINK)  from 
F(I , J)  since  this  additional  flow  along  ARC(I.J)  will 
be  a reversal  of  previous  flow  from  n ■ 0^  to  T^.  The 
next  node  to  consider  is  n *»  T^. 

This  lternatlve  procedure  is  continued  until  n = SOURCE.  At  this  point 
a path  from  the  source  to  the  sink  has  been  retraced  working  backwards 


-26- 


from  the  sink.  The  arcs  on  this  path  that  are  traversed  in  the  forward 
direction  (0^.  -»•  Tj)  as  we  go  from  the  source  to  the  sink  have  their  flows 
increased  by  CAP(SINK)  while  the  arcs  on  this  path  that  are  traversed  in 
the  reverse  direction  (T^  0^)  have  their  flows  decreased  by  CAP(SINK). 

All  labels  are  now  discarded  and  the  labeling  process  (A)  is  started 

over. 

C.  Node  Number  Change.  The  labeling  process  has  resulted  in  non- 
breakthrough. The  following  subsets  of  arcs  are  determined: 

A1  = {(I,J)jOI  labeled,  Tj  unlabeled,  ABAR(I,J)  < 0},  (2.52) 

A2  = {(I,J)|0I  unlabeled,  T labeled,  ABAR(I,J)  > 0).  (2.53) 

We  now  define 

l 


DELTA1  * min[-ABAR(I,J)] , 

(2.54) 

A1 

DELTA2  =*  min[ABAR(I,  J)  ] , 

(2.55) 

A2 

DEL  = rain (DELTAl , DELTA2) . 

(2.56) 

1 

j 


The  node  numbers  XNODE(K)  are  changed  by  subtracting  DEL  from  all  XNODE(K) 
corresponding  to  unlabeled  K.  All  labels  are  discarded  and  the  labeling 
process  (A)  is  started  over. 

2.4.  Flowchart  of  the  Algorithm 


See  Figure  9. 
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FIGURE  9 : Flowchart  of  the  Algorithm 


Read  the  algorithm  parameters: 
NA  = Total  number  of  activities 
NN  = Total  number  of  nodes 
SINK  - Number  of  the  sink  node 
SOURCE  = Number  of  the  source  node 
For  each  activity  I: 


Oj,  Tr  NK(I) 


COST(1 ,1) , . . . .COST [I , NK(I) ] 

TIME (1,1) TIMElI,  NK(I)  ] 

Either  the  entire  project  cost  curve  is 
determined  or  optionally  just  the  optimal 
activity  times  for  one  specified  project 
deadline  time,  LAMBDA. 


I 


Set  all  F(I, J)=0 

Starting  routines  set  all  activity  times  at  their  minimums  and 
determine  the  corresponding  XNODE(K)’s  and  minimum  feasible 
LAMBDA,  LMIN. 

Is  specified  LAMBDA  <.  LMIN? 


YES 


Terminate  the  algorithm. 
No  feasible  activity 
times. 


Set  all  activity  times  at  their  upper 
bounds  and  determine  the  corresponding 
XNODE(K)'s  and  maximum  practical  LAMBDA, 
LMAX. 

Is  specified  LAMBDA  >_  LMAX? 


Terminate  the  algorithm. 
The  optimal  activity  com- 
pletion times  are  their 
upper  bounds. 


Assign  the  SOURCE  the  label  [-,-,-,00], 
all  other  nodes  are  unlabeled. 


Go  toi 


© 


■WPD 
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Find  the  following  subsets: 

A^:  {(I,J)  such  that  Oj  is  labeled,  is  unlabeled,  and  ABAR(I,J)<0} 

A^:  {(I,J)  such  that  0^.  is  unlabeled,  T^.  is  labeled,  and  ABAR(I,J)>0} 

DELTA1  « min[-ABAR(I,J) J 
A1 

DELTA2  - min[ABAR(I, J) ] 

A2 

DEL  = minfDELTAl,  DELTA2 ] . 


Subtract  DEL  from  all  unlabeled  nodes  XNODE(K) . Then  the 
XACT(I)  = min{TIME[I,NK(I)] ,XNODE(TI)-XNODE(0I)} 
are  an  alternative  optimal  solution  for  the  current  LAMBDA  and  also  an 
optimal  solution  for  new  LAMBDA  = current  LAMBDA  - DEL. 

A new  point  on  the  project  cost  curve  has  been  determined. 

Is  new  LAMBDA  specified  LAMBDA? 


Terminate  the  algorithm. 
Desired  solution  found. 


Discard  labels  and  return  to 


©• 
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3.  VERIFICATION  OF  CLAIMS 

The  algorithm  described  In  the  previous  chapter  is  based  on  many 
claims.  The  lemmas  and  theorems  given  below  prove  these  claims  and,  at 
the  same  time,  show  that  the  given  algorithm  does  indeed  yield  for  each 
project  deadline  time  LAMBDA  the  individual  activity  completion  times 

which  minimize  the  total  project  cost. 

- 


i 


' /' 


Summary 


The  initial  primal  and  dual  variables  XACT(I) ,XN0DE(K)  and  F(I,J) 
provide  a feasible  and  optimal  solution  for  the  largest  LAMBDA  of  interest, 
LMAX  (Lemmas  1,  5 and  6).  The  changes  applied  to  these  variables  arise 
from  either  breakthrough  or  nonbreakthrough  and  force  the  variables  to  re- 
main feasible  and  satisfy  the  optimality  properties  throughout  the  algo- 
rithm (Lemmas  2,  3,  4,  5,  6 and  Theorem  1).  The  optimality  properties 
imply  that  complementary  slackness  holds  which,  combined  with  feasibility, 
implies  that  the  solution  is  optimal  for  a given  LAMBDA  (Lemma  7 and 
Theorem  2) . 

The  algorithm  itself  terminates  after  a finite  number  of  applications 
of  the  labeling  procedure  (Theorem  3) . At  the  conclusion  of  the  compu- 
tations a path  from  the  source  to  the  sink  has  been  identified  in  the 
First  Labeling  step  such  that  along  this  path 


TIME(I.l)  + XNODE(0I)  = XN0DE(TI> . 


decrease  in  LAMBDA  would  make  Che  Primal  Problem  infeasible;  l.e.,  the 
project  cannot  be  completed  in  any  shorter  time. 

The  project  cost  function  PCOST (LAMBDA)  is  convex  and  is  linear 
between  the  successively  determined  values  of  LAMBDA  generated  in  the 
computations  (Lemmas  8 and  9).  Given  two  successively  determined  values 
of  LAMBDA,  say  LI  and  L2  * LI  - DEL,  the  optimal  node  times  and  activity 
completion  times  for  any  project  deadline  time  L between  LI  and  L2  are 


XNODE  (K) 

L 

XACT  (I) 
Li 


XNODEl1(K)  if  K labeled  when  LAMBDA=L1 , 

XN0DEL1(K)-(L1-L)  if  K unlabeled  when  LAMBDA=L1 
= min(TIME[I ,NK(I) ] , XNODE  [T_]-XNODEt [0t]} 

Li  X Li  L 


where  the  subscript  LI  implies  LAMBDA  ■ LI  (Theorem  4) . 

One  additional  feature  of  the  algorithm  is  that,  if  the  problem  is 
specified  in  terms  of  integers,  then  the  breakpoints  of  the  project  cost 
curve  PCOST (LAMBDA)  and  the  corresponding  optimal  activity  times  will  all 
be  integers. 


3.2.  Proofs 

Lemma  1:  The  original  set  of  node  integers  XNODE (K)  and  the  zero  flow 
F(I , J)  satisfy  the  optimality  properties.  Furthermore,  this  F(I,J)  mini- 
mizes (2.36)  implying  an  optimal  solution  for  LAMBDA  = LMAX. 

Proof : In  a starting  routine  the  activity  times  XACT(I)  are  set  to 
their  largest  feasible  (cheapest)  values.  Then  the  node  times  XNODE(K) 
are  set  to  their  corresponding  smallest  feasible  values.  This  implies 


that 
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TIMEfl,  NK(I) ] <_  XNODE (T^)  - XNODE(Ot) 


or  equivalently 

TIME[I,  NK(I)  ] + XNODE (0^)  - XNODE^)  <_  0. 

Thus  all  ABAR(I,  J)  _<  0.  Finally,  since  all  ABAR(I,J)  < 0 and  F(I,J)  * 0, 
the  optimality  properties  are  satisfied. 

The  dual  objective  function  is 

LAMBDA  • V - l TIME(I,NK(I)  + 1 - J)  • F(I,J) 

I,J 

“ -t  I TIME(I,NK(I)  + 1 - J)  • F(I , J)  - LAMBDA  • V] 

I.J 

- -{  I TIME(I,  NK(I)  + 1 - J)  • F(I, J)  + [XNODE (SOURCE) 

I.J 

- XNODE(SINK) ] • V) 


- -{  l TIME(IfNK(I)  + 1 - J)  • F(I , J)  + l [ XNODE (0.) 

I.J  I.J  1 

- XNODE (TT)1  * F(I , J) } 

- -l  I ABAR(I.J)  * F(I.J)]. 

I.J 

Thus,  since  all  ABAR(I.J)  £ 0,  F(I,J)  = 0 is  optimal.  QED. 

Lemma  2:  If  breakthrough  occurs,  the  old  node  numbers  and  the  new  flow 
satisfy  the  optimality  properties. 

Proof:  The  node  numbers  XNODF.(K)  do  not  change.  The  new  flows  arc 
obtained  by  adding  the  positive  number  CAP  (SINK)  to  all  F(I,J)  correspond- 
ing to  forward  arcs  of  the  path  from  the  source  to  the  sink,  and  subtract- 
ing CAP(SINK)  from  all  F(I,J)  corresponding  to  reverse  arcs  of  the  path. 


| ■■ 
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Flow  changes  occur  only  In  arcs  for  which  ABAR(I,J)  = 0.  No  restriction 
is  imposed  on  the  F(I,J)'s  in  the  optimality  properties  when  ABAR(I,J)  * 0. 
Thus,  the  old  XN0PE(K)'s  and  the  new  F(I,J)'s  still  satisfy  the  optimality 
properties.  QED. 

Lemma  3:  If  nonbreakthrough  occurs,  the  node  number  change,  DEL,  is  a 
well-defined  positive  number. 

Proof:  For  DEL  to  be  well-defined,  at  least  one  of  the  sets  of  arcs 
A^,  A^  (as  defined  in  equations  (2.52)  and  (2.53))  is  non-empty. 

Suppose  A^  were  empty.  Since  there  is  a path  from  the  source  to  the 
sink  in  the  project  network,  and  since  the  source  is  labeled  and  the  sink 
is  unlabeled,  there  must  be  a set  of  arcs  (ARC(I,J),  J = 1,  ...,  *TK(I)} 
in  the  enlarged  network  with  0^  labeled  and  unlabeled.  The  definition 
of  A^  implies  that  if  A^  is  empty,  then  ABAR(I,J)  >_  0 for  this  set  of 
arcs.  From  labeling  rules  (2.46)  and  (2.47),  if  ABAR(I,J)  = 0 then  F(I,J) 
cannot  be  less  than  (C[I,  NK(I)  - J]  - C[I,  NK(I)  - J +1]},  otherwise 
Tj  would  have  been  labeled  from  0^.  From  (2.43),  if  ABAR(I,J)  > 0,  this 
implies  that  F(I,J)  = C[I,  NK(I)  - J]  - C[I,  NK(I)  - J + 1].  Hence'we 
have  F[I,  NK(I)]  = <=°.  But  this  F[I,  NK(I)]  is  part  of  the  actual  flow  through 
the  network  and,  if  it  equals  infinity,  the  first  labeling  process  would 
have  terminated  the  algorithm.  Since  this  has  not  happened,  there  are  no 
infinite  flows  and  A^  is  non-empty. 

By  definition,  DEL  is  always  positive.  QED. 

Lemma  4:  If  nonbreakthrough  occurs,  for  any  DEL'  satisfying  0 <_  DEL'  _<  DEL, 
the  new  node  numbers 


■ 1 
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XNODE' (K) 


fxNODE(K) 


l> 


XNODE(K)  - DEL' 


for  K labeled, 
for  K unlabeled, 


. I 


and  the  old  flow  F(I,J)  still  satisfy  the  optimality  properties. 

Proof:  The  new  AMR'(I.J)  = TIME(I,NK(I)  + 1 - J)  + XNODE' (0.^  - XNODE' (Tj), 
(i)  Suppose  ABAR'(I,J)  < 0.  Then  F(I,J)  = 0 because  of  the 

following: 

(a)  If  A3AR(I, J)  < 0,  then  F(I,J)  * 0 by  (2.42). 

(b)  If  ASAR(I,J)  = 0,  then 

TIME  (I  ,NK(I)  + 1 - J)  + XNODE  (0^.)  - XNODE^)  * 0, 

or  equivalently 

TIME  (I  ,NK(I)  + 1 - J)  = - XNODECOj)  + XNODECTj.)  ; 

so  that 

ABAR' (I , J)  = TIME(I,NK(I)  + 1 - J)  + XNODE' (O^.)  - XNODE' (T^)  < 0, 

implies 

- XNODE (Oj)  + XNODE (Tj.)  + XNODE'  (0^.)  - XNODE'  (Tj.)  < 0, 

and  finally 

XNODE' (0j)  - XNODE (Oj)  < XNODE' (T^)  - XNODECTj.); 
but  this  can  happen  only  when  0^  is  unlabeled  and  T^  is  labeled.  Hence, 
if  ABAR(I,J)  = 0,  then  by  labeling  rules  (2.49)  and  (2.50),  F(I,J)  = 0, 
otherwise  0^  would  be  labeled  from  T^.. 

(c)  If  ABAR(I , J)  > 0,  then 

TIME(I ,NK(I)  + 1 - J)  + XN0DE(0].)  - XNODE(Tj)  > 0, 

or  equivalently 

TIME(I ,NK(I)  + 1 - J)  > - XNODE(0].)  + XNODE (Tj.)  ; 

so  that 


: 


-35- 


lmplies 


and  finally 


ABAR' (I, J)  = TIME(I,NK(I)  + 1 - J)  + XNODE' (Oj.)  - XNODE' (Tj)  < 0, 


TIME(I,NK(I)  + 1 - J)  < XNODE'  (TJ)  - XNODE' (Oj.), 


XNODE ' (T^)  - XNODE' (Oj)  > - XNODECOj)  + XNODE (Tj. ) , 


XNODE'  (T^.)  - XNODE(T^)  > XNODE’  (Oj)  - XNODECO^.) 


Again,  this  can  happen  only  when  Oj  is  unlabeled  and  T is  labeled.  But 
then  the  arc  ARC(I,J)  is  in  and  DEL  <_  ABAR(I,J).  This  would  imply 


ABAR' (I, J)  = TIME(I,NK(I)  + 1 - J)  + XNODE(Ox)  - DEL  - XNODE (T^) 
= ABAR ( I , J)  - DEL 


which  contradicts  the  assumption  ABAR' (I, J)  < 0.  Hence  this  case  cannot 


(ii)  Suppose  ABAR' (I, J)  = 0.  There  are  no  restrictions  on 
F(I,J)  so  the  optimality  properties  still  hold. 

(iii)  Suppose  ABAR' (I, J)  > 0.  Then 


F(I , J)  » Cfl,  NK(I)  - J]  - C[I,  NK(I)  - J + 1] 


because  of  the  following: 


(a)  If  ABAR(I,J)  > 0,  F(I,J)  - C[l,  NK(I)  - Jj  - C[I,  NK(I)  - J + 1] 


by  (2.43). 


(b)  If  ABAR(I, J)  - 0,  then 


ABAR(I , J)  < ABAR' (I, J) 


or  equivalently 


TIME(I,NK(I)  + 1 - J)  + XNODEfOj)  - XNODE (T^) 

< TIME  (I  ,NK(I)  + 1 - J)  + XNODE'COj)  - XNODE '(Tj.); 

so  that 

XNODECOj)  - XNODE' (Oj)  < XNODE (T^)  - XNODE' (Tj). 

This  can  happen  only  if  0^  is  labeled  and  T^  is  unlabeled.  Hence,  by 
labeling  rules  (2.46)  and  (2.47) 

F(I, J)  = C[I,  NK(I)  - J]  - C[I,  NK(I)  - J + L], 

otherwise  T^.  would  be  labeled  from  0j. 

(c)  If  ABAR(I, J)  < 0,  then 

ABAR(I,J)  < ABAR'(I.J), 

and  again 

XNODE(0I)  - XNODE'  (Oj)  < XNODE  (T^.)  - XNODE'  (T^)  . 

This  can  happen  only  if  Oj  is  labeled  and  T^  is  unlabeled.  But  then  the 
arc  ARC(I , J)  is  in  A^  and  DEL  <_  - ABAR(I,J)  which  would  imply  that 

ABAR'(I.J)  = TIME (I ,NK(I)  + 1 - J)  + XNODE' (.0^)  - XNODE' (T^) 

= TIME(I ,NK(I)  + 1 - J)  + XNODE (Oj)  - XNODE^)  + DEL 

= ABAR(I , J)  + DEL 

£ 0. 

This  contradicts  the  assumption  ABAR'(I,J)  > 0.  Hence  this  case  cannot 
occur. 

Cases  (i)  - (iii)  together  imply  that  the  new  node  numbers  and  the 

old  flow  still  satisfy  the  optimality  properties.  QED. 
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Theorem  1:  The  optimality  properties  (2.42)  and  (2.43)  are  maintained 
throughout  the  algorithm. 

Proof : From  I.emma  1,  we  see  that  the  initial  node  numbers  XNODE(K)'s 
and  the  zero  flow  provide  an  optimal  solution  for  LAMBDA  = LMAX.  If  break- 

7 

through  occurs,  we  see  that  the  new  F(I,J)'s  are  still  optimal  (Lemma  2). 

If  nonbreakthrough  occurs,  we  have  a well-defined  positive  number  DEL  with 
which  to  update  the  XNODE(K)'s  (Lemma  3)  and,  from  Lemma  4,  these  updated 
values  satisfy  the  optimality  properties.  QED. 

Lemma  5:  The  starting  values  of  the  XNODE(K)'s  and  XACT(I)'s  are  feasible 
and  remain  feasible  throughout  the  algorithm. 

Proof ; The  starting  values  are  found  by  an  algorithm  that  sets  the 
XACT (I) ' s to  their  largest  feasible  times,  TIME[I,  NK(I)].  Correspondingly 
the  XACT(I,M)'s  are  set  equal  to  their  upper  bounds  and  hence  (2.15)  and 
(2.16)  are  satisfied.  Then  the  algorithm  sets  XNODE(K)  equal  to  the 
length  of  the  longest  path  from  the  source  to  node  K,  which  implies  that 
(2.13)  is  satisfied.  We  also  define  XNODE(SOURCE)  = 0 and  LAMBDA  = 
XNODE(SINK);  hence  (2.14  ) is  satisfied  and  the  initial  values  are  feasible. 

If  breakthrough  occurs,  the  XNODE(K)'s  and  XACT (I) 's  are  not  changed 
and  hence  remain  feasible. 

If  nonbreakthrough  occurs,  the  labeled  XNODE(K)'s  are  unchanged,  and 
the  unlabeled  XNODE(K)'s  are  updated  by  subtracting  DEL,  determined  by  (2.56) 


new  LAMBDA  = XNODE(SINK)  - DEL 


so  (2.14)  is  satisfied. 

(i)  Suppose  both  Oj  and  Tj  are  labeled  for  activity  I.  Then 

neither  these  nodes  nor  XACT(I)  are  updated  and  hence  XACT(I)  remains 


Suppose  both  0 and  T are  unlabeled  for  activity  I.  Then 


new  XNODE (0)  - old  XNODE (0T)  - DEL 


new  XNODE (T_)  - old  XNODE(Tt)  - DEL,  and 


new  XACT(I)  - min{TIME[I,  NK(I)],  old  XNODE(Tx)  - DEL 


old  XNODE (0T)  + DEL} 


old  XACT(I) 


£ old  XNODE (T_)  - old  XNODE (0T) 


new  XNODE (T_)  - new  XNODE (0T) 


or  equivalently 


new  XACT(l)  + new  XNODE(0_)  - new  XNODE(T_)  £ 0 


so  that  (2.13)  Is  satisfied.  Since  XACT(I)  has  not  changed,  (2.15)  and 


(2.16)  are  still  satisfied.  Therefore,  in  this  case,  feasibility  is 


maintained 


(iii)  Suppose  0 is  labeled  and  Tt  is  unlabeled  for  activity  I 


new  XNODE (Tt)  = old  XNODE (Tt)  - DEL 


new  XNODE (0)  ■ old  XNODE (0T ) , and 


new  XACT(I)  - min{TIME[I,  NK(I)J,  old  XN0DE(Tt)  - DEL 


old  XNODE (0_)} 
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so  that  (2.13)  and  (2.15)  are  satisfied.  The  lower  bound  constraint, 
(2.16),  is  also  satisfied  because  of  the  following: 

(a)  Suppose  ABAR[I,NK(I) ] < 0.  Then  since  0^.  is  labeled  and  T^  is 
unlabeled,  the  definition  of  DEL  implies  that 


XN0DE(TI)  - XNODE(0].)  - TIME (1,1)  _>  DEL 


and  hence 


XN0DE(Tt)  - XNODE(0].)  - DEL  >_TIME(I,1) 


which  implies  that  XACT(I)  >_TIME(I,1). 


(b)  Now  ABAR[I,  NK(I)  ] = 0 cannot  occur,  since  T^.  would  have  been 
labeled  from  0^. 

(c)  Also  ABAR[I,  NK(I)]  > 0 cannot  happen  since  this  would  imply  that 


old  XNODE(0].)  + TIME(1, 1)  _>  old  XNODEOTj) 


which  contradicts  the  feasibility  of  the  previous  node  times. 

(iv)  Suppose  0j  is  unlabeled  and  T^  is  labeled  for  activity  I. 


new  XNODE(0I)  = old  XNODE(Or)  - DEL  , 


new  XNODEd'j)  = old  XNODEdj),  and 


new  XACT(I)  - min{TIME[I,  NK(I)],  old  XNODE^) 


- old  XNODE(0].)  + DEL}; 


so  that  (2.13)  and  (2.15)  are  satisfied.  Since 


TIME(I,1)  ^ old  XACT(I)  £ new  XACT(I), 


the  lower  bound  constraint,  (2.16),  is  trivially  satisfied.  QED. 

Lemma  6;  The  starting  values  of  the  F(I,J)'s  and  V are  feasible  and  re- 
main feasible  throughout  the  algorithm. 

Proof:  Initially,  the  values  of  the  F(I,J)'s  and  V are  set  to  zero. 
Conservation  of  flow,  (2.24),  is  trivially  satisfied  since  the  flow  going 
into  each  node  is  equal  to  zero  which  is  also  equal  to  the  flow  going  out 
of  each  node,  l.e. 

I [If<I,J)]=0«  l [If(I,J)]. 

IeO^K  J IeTj-K  J 

Since  all  F(I,J)  are  set  equal  to  zero,  they  satisfy  their  upper  and 
lower  bounds.  Hence,  the  starting  values  are  feasible. 

If  nonbreakthrough  occurs,  the  values  for  F(I,J)  do  not  change, 
hence  remain  feasible. 

If  breakthrough  occurs,  the  F(I,J)  along  the  path  from  the  source  to 
the  sink  are  updated  by  a positive  number  CAP (SINK)  determined  by  (2.48) 
or  (2.51);  all  other  flows  remain  unchanged.  Suppose  activity  I is  an 
arc  along  the  path  from  the  source  to  the  sink.  Then  either  T^  is  labeled 
from  Oj  or  0^  is  labeled  from  Tj. 

(i)  In  the  former  case,  F(I,J)  < C[I,  NK(I)  - J]  - C[I,NK(I) 

-J  + 1]  by  labeling  rules  (2.46)  and  (2.47),  and  CAP(I)  is  given  by 
(2.48).  This  CAP(I)  is  the  minimum  of  the  previous  CAP  and  C[I,  NK(I) 

- J]  - C[I,  NK(I)  - J + 1]  - F(I , J)  > 0.  Now  CAP (SINK)  < CAP(I) 
and  new  F(I,J)  - old  F(I,J)  + CAP(SINK).  Conservation  of  flow  is  satis- 
fied since  the  same  value  is  added  to  or  subtracted  from  all  activities 
along  this  path  and  V,  the  total  flow,  is  increased  by  CAP(SINK).  Also, 
(2.34)  is  satisfied  since 
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, 

i> 


I 


I 

0 < old  F(I,J)  + CAP (SINK)  < old  F(I,J)  + CAP(I) 

< C(I,  NK(I)  - J]  - C{I,  NK(I)  - J + U. 

Hence  the  new  F(I,J)'s  are  feasible. 

(11)  In  the  latter  case.  F(I,J)  > 0 by  labeling  rules  (2.49) 

and  (2.50),  and  CAP(I)  Is  given  by  (2.51);  i.e.,  the  minimum  of  the  pre 
vlous  CAP(K)  and  F(I,J).  Conservation  of  flow  is  again  satisfied.  The 
following  also  shows  that  (2.34)  is  satisfied:  Now 

old  F(I, J)  < C[I,  NK(I)  - J]  - C[I,  NK(I)  - J + 1],  and 

new  F(I, J)  - old  F(I,J)  - CAP(SINK). 

Since  CAP (SINK)  < CAP(I)  <_  old  F(I,J),  this  implies  that 

0 £ new  F(I,J)  < C[I,  NK(I)  - J]  - C^T,  NK(I)  - J + 1], 

Hence,  F(I,J)  remains  feasible  for  this  case  as  well  and,  therefore,  re- 
mains feasible  throughout  the  algorithm.  QED. 

Lemma  7:  The  optimality  properties  (2.42)  and  (2.43)  imply  that  comple- 
mentary slackness  holds  between  the  primal  and  the  dual  problems. 

Proof:  We  will  use  the  original  pair  of  primal  and  dual  problems 
((2.13)  - (2.19)  and  (2.22)  - (2.25)  respectively)  along  with  the  defini- 
tions of  G(I,M),  H(I,M)  and  F(I,J)  to  show  that  the  complementary  slack- 
ness conditions  are  satisfied;  l.e., 

(I)  XACT(I)  + XNODE(0I)  - XN0DE(T].)  < 0 
implies  that  F(I)  - 0; 

(II)  XACT(I,M)  < U(I,M)  implies  that  G(I,M)  - 0;  and 

(ill)  XACT(I,M)  > L(I,M)  implies  that  H(I,M)  - 0. 
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; 

(i)  If  XACT(I)  + XNODE(Oj)  - XNODE(Tj)  < 0,  then 

XACT(I)  < XNODE(Tx)  - XNODECOj). 

Since 

XACT(I)  - min{TIME[I,  NK(I)],  XNODE(Tj)  - XNODECOj)}, 
this  implies  that 


XACT(I)  - TIME [I » NK(I)]. 


Hence , 

TIME  [I,  NK(I)  ] + XNODECOj.)  - XNODEdj)  < 0. 


Since 


TIME (I , 1)  < TIME(I,  2)  < ...  < TIME [I,  NK(I)], 


it  follows  that 

TIME(I,  NK(I)  + 1 - J)  + XNODE(0].)  - XNODE^)  < 0 

for  J - 1,  2,  . ..,  NK(I).  From  optimality  property  (2.42),  F(I,J)  - 0 for 

J - 1,  2,  ....  NK(I)  , and  finally  F(I)  = £f(I,J)  = 0. 

J 


(Remark:  Since 

ABAR(I,J)  =>  TIME  [I,  NK(I)  + 1 - J]  + XNODE^)  - XNODE(Tj) 

and 


TIME(I.l)  < TIME(I,2)  < 


< TIME[I,  NK(I) ] , 


AJ . 


• • • 
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it  follows  that 

ABAR(I,1)  ^ABARfI, 2)  >_  ...  _>  ABARfI,  NK(I)]. 

Now  the  TIME(I,J)'8  will  be  strictly  Increasing  and  the  ABAR(I,J)'8  strictly 
decreasing  unless  there  is  only  one  possible  value  for  XACT(I)  in  which 
case  the  upper  and  lower  bounds  for  F(I)  and  the  F(I,J)'s  are  0. 

Therefore  in  the  Second  Labeling  part  (i) , page  24,  there  can  only 
be  one  J such  that 

ABARfI, J)  - 0 


and 


F(I , J)  < C[I,  NK(I)  - J]  -C[I,  NK(I)  - J + 1]. 


For  this  J ’ 

0 > ABARfI, J +1)  > ...  > ABARfI,  NKfl) ] , 
so  that  by  optimality  property  f2.42)  ' 

Ffl,  J + 1)  = ...  - F[I,  NKfl)]  = 0. 

Also,  for  this  J 

I 

ABARfI,  1)  > ...  > ABARfI, J - 1)  > 0, 

so  that  by  optimality  property  (2.43)  Ffl,  1)  ....  Ffl,  J - 1)  are  all  at 
their  upper  bounds.  Thus,  when  Ffl)  is  increased,  it  is  the  F(I,J)  with 
the  smallest  index  J such  that  F(I,J)  is  less  than  its  upper  bound  which 


is  increased. 


Similarly  the  Second  Labeling  part  (ii)  and  the  optimality  properties 
imply  that  when  F(I)  is  decreased  it  is  the  F(I,J)  with  the  largest  index 
J such  that  F(I,J)  > 0 which  is  decreased.  Therefore,  if  F(I,J)  is  posi- 
tive,then  F(I,1) , . . ,F(I, J - 1)  are  all  at  their  upper  bounds;  and,  if 
F(I,J)  * 0,  then  F(I,  J + 1) F[I,  NK(I)]  also  equal  0.  These  natur- 

al properties  of  the  F(I,J)’s  are  used  in  parts  (ii)  and  (iii)  below.) 


Show  that  XACT (I ,M)  < U(I,M)  implies  G(I,M)  - 0 where,  as 


in  (2.15)  and  (2.28), 


:,M)  =| 


TIME(I,2) 

TIME(I,M  + 1)  - TIME(I,M) 


M = 1 


M - 2,  ....  NK(I)  - 1, 


G(I,M)  = max{ 0 , C(I,M)  - F(I)}, 


XACT (I, M) 


•min[U(I,M),  XACT (I)] 


M = 1 


min{U(I,M),  max[0,  XACT(I)  - TIME(I,M)]} 


M = 2,  ....  NK(I)  - 1. 


If  XACT ( I , M)  < U(I,M),  then 


XACT(I,M) 


XACT (I) 

max[0,  XACT ( I ) - TIME(I,M)] 


M = 1 


M = 2,  ...»  NK(I)-1. 


Case  I : M - 1. 


Since 


TIME(I,  2)  = U(l,  1)  > XACT (I , 1)  - XACT(I), 


It  follows  that 


TIME(I,  2)  > XACTCI)  - min{TIME[I,  NK(I)],  XNODE(T_)  - XN0DE(0_)} 


TIME (I , 2)  < TIME [I,  NK(I) ] 


this  Implies  that 


TIME(I,  2)  > XNODE(T_)  - XNODE(0_) 


ABAR [ I , NK(I)  - 1]  - TIME(I,2)  + XNODE(0T)  - XNODE(Tx)  > 0 


Therefore  F(I,J)  *=  C[NK(I)  - J]  - C[NK(I)  - J + 1] 


Hence 


NK(I)-1 


F[I,  NK(I) ] + C[I,  NK(I)  - 1]  - C[I,  NK(I) ] 


F[I,  NK(I) ] + C(I,  1)  - C [I , NK(I) ] 


Since  ClI,  NK(I)]  = 0,  F(I)  >_  C(I,  1).  Therefore 


0. 


Now  U ( I , M)  > XACT(I.M)  - max[0,  XACT(I)  - TIME(I,M)]  implies 


U(I,M)  > XACT(I)  - TIME(I,M) 


U(I,M)  > min{TIME[I,  NK(I)],  XNODE(T  ) - XNODE(Ot)}  - TIME(I.M) 


U(I,M)  + TIME(I.M)  > min{TIME[I.  NK(I)],  XNODE(T  ) - XNODE(0  )} 


Since  U(I,M)  * TIME(I,M  + 1)  - TIME(I.M) 


U(I ,M)  + TIME(I,M)  = TIME(I,  M + 1) 


TIME( 


1)  > min{TIME[I,  NK(I)],  XNODE(T_)  - XNODE(0_)} 


Since  TIME(I,M  + ^ TIME [I,  NK(I)],  this  implies 


TIME(I,  M + 1)  > XNODE(T  ) - XN0DE(0_) 


TIME(I,M  + 1)  + XNODE(0_)  - XNODE(T_)  > 0 


F[I,  NK(I)  - M - 1]  are  also  at  their  upper  bounds 


C[I * NK(I)  - 1]  - C[I,  NK(I) ] 


C(I,  M)  - C[I,NKCD] 


NK(I)  - 1 


(iii)  Show  that  XACT(I,M)  > L(I,M)  implies  that  H(I,M)  = 0 where,  as 


TIME(I,  1) 


min[U(I,M) , XACT(I))  M - 1 

XACT(I,M)  - < 

min{U(I,M) , max[0,XACT(I)  - TIME(I,M)]} 

M * 2,  ....  NK(I)-1. 

Case  I;  M ■ 1. 

If  XACT(I,I)  > L(I ,1) , then 

TIME(I.l)  - L(I,1)  < XACT(I,1)  - min[U(I,l),  XACT(I)]. 

Since  0(1,1)  - TIMS(I,2), 

TIME(I.l)  < aln [TIME (I ,2),  XACT(I)J 

and 

TIME(I , 1)  < XACT(I). 

Thus 

TIME(I , 1)  < XACT(I)  = min{TIME[I , NK(I)],  XNODE(T^)  - XNODECOj)} 

and 

TIME (1,1)  < XNODE(T  ) - XNODE(Ot). 

1 1.  » ' 

Therefore 


ABAR[I,  NK(I) ] - TIME(I.l)  + XNODE(0_)  - XN0DE(T_)  < 0. 


+ C [I » NK(I)  - 2]  - C[I,  NK(I)  - 1] 


H(I, 1)  = max[0,F(I)  - C(I,1)] 


0 * L(I,M)  < XACT (I ,M)  »*  min{U(I,M) .max [0, XACT (I)  - TIME(I,M)]} 


0 < min{U(I,M) , max[Q,  XACT(I)  - TIME(I,M)]} 


0 < max[0,  XACT (I)  - TIME(I,M)],  and 


0 < XACT (I)  - TIME (I ,M) 


This  implies  that 


TIME(I,M)  < XACT (I)  =■  min {TIME II,  NK(I)],  XNODE(T  ) - XNODE(0  )} 


TIME(I.M)  < XNODE (T  ) - XNODE(Ox) 


ABAR [ I , NK(I)  - M + 1]  - TIME(I,M)  + XNODE(0  ) - XNODE (T)  < 0 


C(I,M)  - C[I,  NK(I) ] 


Theorem  2:  Since  the  XNODE(K)'s,  XACT(I,M)'s,  V,  AND  F(I) 'a  are  feasible  and 


complementary  slackness  holds,  they  are  optimal 


Proof:  The  primal  problem  (2.13)  - (2.19)  is  in  the  form 


subject  to 


where  the  x vector  contains  the  XNODE(K)'s  and  XACT(I,M)'s.  The  dual  problem 


where  the  w vector  contains  V and  the  F(I)'s 


holds  for  any  feasible  x and  w.  When  Ax  < b is  rewritten  in  the  form 


?! 


/ 


where  xg  is  a vector  of  slack  variables,  complementary  slackness  means 


w x * 0. 
s 


Therefore,  since  for  any  feasible  x and  w 


T T T 

wAx  + wx  **  w b 
s 


T T T 

cx  + wx  = b w„ 
s 


complementary  slackness  Implies 


T . T 
c x = b w 


and  hence  that  both  x and  w are  optimal.  QED. 

Theorem  3:  The  algorithm  terminates  after  finitely  many  applications  of  the 
labeling  procedure. 

Proof:  In  order  that  the  algorithm  fail  to  terminate,  an  infinite  sequence 
of  breakthroughs  and  nonbreakthroughs  would  have  to  occur. 

I 

Since  the  flow  change  following  a breakthrough  has  a positive  minimum,  an 
infinite  number  of  breakthroughs  would  produce  flows  having  arbitrarily  large 
values  V.  However,  when  a sufficiently  large  value  V is  reached,  there  will  be 
a path  from  the  source  to  the  sink  with  F[I,  NK(I) ] > 0 all  along  this  path. 

Since  ABARfl , NK(I) ] <_  0 throughout  the  computations,  we  would  have  ABAR[I ,NK(I) ] - 0 
for  arcs  on  this  path.  But  then  the  first  labeling  procedure  would  terminate. 
Therefore,  there  can  only  be  a finite  number  of  breakthroughs. 

Following  a nonbreakthrough,  all  nodes  previously  labeled  can  again  be 
labeled.  (This  follows  from  the  fact  that  for  labeled  0^  and  T^ , the  new  ABAR(I,J) 
is  equal  to  the  old  ABAR(I,J)).  In  addition,  at  least  one  more  node  can  be 


labeled  (the  node(s)  corresponding  to  the  arc(s)  In  A 


DEL).  Eventually,  the  number  of  nodes  that  can  be-  labeled  will  reach  the  total 


number  of  nodes  Implying  that  the  sink  can  be  labeled  and  the  occurrence  of  a 
breakthrough.  Therefore,  Infinitely  many  successive  nonbreakthroughs  cannot  occur 
Hence,  there  can  only  be  a finite  number  of  applications  of  the  labeling 
procedure.  QED. 


Definition:  A function  P(X)  Is  said  to  be  convex  over  some  interval  in  X,  if 


for  any  two  points  XI,  X2  in  the  Interval  and  for  all  a,  0 < a < 1 


Lemma  8:  PCOST (LAMBDA)  Is  convex  for  LMAX  > LAMBDA  > LMIN,  where 


LMAX  s the  longest  (cheapest)  time  to  complete  the  project 


LMIN  = the  shortest  time  to  complete  the  project 


Proof:  Let  LI  > L2  both  be  in  the  interval  [LMIN,  LMAX].  Let 


for  some  a in  [0,  1].  Also  let  XACT  (I) , XNODE. (0_)  XNODE. (T  ) , XACT.(I) 


XNODE»(0  ),  and  XN0DE„(TT)  represent  optimal  solutions  to  the  problems  corres 


ponding  to  LAMBDA  ■ LI  and  LAMBDA  - L2  respectively.  We  first  want  to  show  that 


[otXNODE  (K)  + (1  - a)  XNODE' (K)  ] and  [a-XACT  (I,M)  + (1  - a)  XACT  (I,M)]  are 


feasible  when  LAMBDA  ■ L.  This  result  follows  easily  since  the  constraints 
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Since  XNODE^Oj)  + JXACT^I.M) 
M 


- XNODE1(TI)  £ 0 


XNODE2(Oi)  + JxACT2(I,H)  - XNODE2(Ti)  £ 0, 
M 


it  follows  that 


[aXN0DE,(0_)  + (1  - a)  XN0DE1(0T)]  + J[aXACT  (I,M)  + (1  - a)  XACT-U.M)] 
2 1 A 1 M 

- [aXNODE^Tj)  + (1  - a)  XNODE^Tj)] 


and  the  constraints  (2.13)  are  satisfied. 

(ii)  Now  XNODE^SINK)  < LI  and  XNODE2(SINK)  £ L2;  so  that 

(1  - o)XNODE1(SINK)  + aXNODE2 (SINK)  £ (1  - a)Ll  + oL2  - L 


and  constraint  (2.14)  is  satisfied. 


(iii)  Also,  L(I,M)  £ XACT1(I,M)  £ U(I,M)  and  L(I,M)  £ XACT2(I,M)  £ U(I,M) 


implies 


L(I,M)  £ (1  - a)XACT1(I,M)  + aXACT2(I,M) 


< U(1,M) 


and  hence  constraints  (2.15)  and  (2.16)  are  satisfied. 
Recall  that 

PCOST (LAMBDA)  - KK  - l [C(I,M)XACT(I,M) ] 

I.M 

where 
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KK  - ][ [COST (I ,1)  + C (1,1) TIME (1,1) ] . 


Hence , 


ctPC0ST(L2)  + (1  - a) PCOST (LI) 

- a(KK  - l [C(I,M)XACT2(I,M)]}  + (1  - a) {KK  - l [C(I ,M)XACT  (I ,M) ] } 

I,M  I,M 

- aKK  + (1  - o)KK  - a £ C(I,M)XACT.(I,M)  - (1  - a)  l C(I,M)XACT  (I,M) 

I,M  1 I ,M  1 

- KK  - l C(I,M)[oXACT  (I,M)}-  I C(I,M)f(l  - a)XACT  (I,M)] 

I ,M  I,M 

- KK  - l C(I,M)[aXACT  (I,M)  + (1  - a)XACT  (I,M)]. 

I ,M  1 


Furthermore  aPC0ST(L2)  + (1  - a)PCOST(Ll)  is  the  objective  function  value  corres- 
ponding to  [aXNODE2(K)  + (1  - a)XNODE1(K)]  and  [aXACTjd.M)  + (1  - oOXACT^I  ,M)  ] 
which  we  have  just  shown  are  feasible.  Therefore,  since  we  are  minimizing 
PCOST (LAMBDA) , 

PCOST(L)  PCOST [aL2  + (1  - a) LI] 

£ aPC0ST(L2)  + (1  - a)PCOST(Ll) , 
and  PCOST (LAMBDA)  is  convex. 

Lemma  9:  The  project  cost  function,  PCOST (LAMBDA) , is  piecewise  linear. 

Proof:  Let  LI  > L2  * LI  - DEL  be  two  successively  determined  LAMBDA'S 
where  DEL  is  determined  by  (2.56).  (Of  course,  LI  could  be  the  initial  value 
of  LAMBDA.)  Suppose  LI  _>  LAMBDA  _>  L2  and  that  F(I,J)'s  and  V were  the  flows 
when  LAMBDA  was  changed  from  L2  to  LI.  Recall  that 

PCOST (LAMBDA)  - PCOST [XNODE( SINK) ] - [[KK(I)  - £c(I ,M)XACT(I ,M) ] 

I M 
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fc'H 
r ; 


which  Is  the  primal  objective  function.  Since  the  primal  and  dual  objective 
functions  are  equal  under  optimality,  we  have,  for  all  LAMBDA  with  LI  _>  LAMBDA 
_>  L2, 

PCOST (LAMBDA)  « Z - LAMBDA* V + l F(I , J)TIMElI , NK(I)  - J + 1] 

IJ 

where  Z is  a constant.  Therefore 


PCOST (LAMBDA)  - PCOST (Ll) 

* Z - LAMBDA *V  + [ F(I, J)TIME[I,  NK(I)  - J + 1] 

IJ 

- Z + Ll* V - £F(I,J)TIME[I,  NK(I)  - J + 1] 
IJ 

- - LAMBDA *V  + Ll* V 


- (Ll  - LAMBDA) *V 

for  Ll  _>  LAMBDA  ^ L2,  so  that  PCOST  (LAMBDA)  is  linear  on  the  given  interval. 
QED. 

Theorem  4;  If  Ll  > L2  = Ll  - DEL  are  two  successively  determined  values  of 
LAMBDA  where  DEL  is  determined  by  (2.56),  then  for  any  value  of  L such  that 
Ll  > L ^ L2  the  optimal  values  of  the  XNODE(K)'s  and  XACT(I)'s  for  that  L are 
given  by 


XNODE  (K) 

L 


XNODEL1(K) 


XN0DEl1(K)-(L1-L) 


if  K is  labeled  when 
LAMBDA  - Ll, 
if  K is  unlabeled  whti 
LAMBDA  - Ll, 


XACTl(I)  - min(TIME[I,NK(I) ] , XNODE^Tj)  - XNODE^Oj)} 


i 


where  the  subscripts  L and  LI  Imply  LAMBDA  ■ L and  LI  respectively. 

Proof:  Since  Lemma  1 states  that  we  begin  with  an  optimal  solution  when 
LAMBDA  - LMAX,  we  can  without  loss  of  generality  assume  that  we  have  found  op- 
timal solutions  for  all  LAMBDA  values  produced  by  the  nonbreakthrough  procedure 
up  to  LAMBDA  - LI.  We  will  now  show  that  the  above  XNODE(K) 's  and  XACT(I)*s 
are  optimal  for  all  LAMBDA  between  LI  and  L2  including  L2.  The  terms  "labeled" 
and  "unlabeled"  below  refer  to  "labeled  when  LAMBDA  - LI"  and  "unlabeled  when 
LAMBDA  = LI"  respectively. 

We  first  want  to  show  that  for  Ll  > L £ L2  the  XNODE  (K) ' s and  XACTt(I)’s 

L L 

are  feasible.  Since  the  definition  of  XACTL(I)  implies  that 
XACTl(I)  + XNODE^Oj.)  - XNODE^Tj.)  £ 0 


XACTl(I)  £ TIMEfI,  NK(I) ] , 

(2.13)  and  (2.15)  are  satisfied.  Therefore,  the  only  aspect  of  feasibility 
left  to  show  is  (2.16),  i.e. 


TIME(I.l)  £ XACT  (I) 


or  equivalently 


TIME(I,1)  £ XNODEL(Tx)  - XNODE^Oj). 


(i)  Suppose  Oj  and  T^  are  both  labeled  for  a specific  activity  I.  Then 
XNODEL(TI)  - XNODEl(Oi)  - XNODELl(TI)  - XN0DELl(0T)  £ TIME(I.l) 


since  the  solution  at  Ll  is  feasible. 
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(il)  Suppose  Oj.  is  labeled  and  is  unlabeled.  Then 

XNODE^Tj)  - XNODE^Oj)  - XNODE^Tj)  - (LI  - L)  - XNODF^Oj). 

The  definition  of  DEL  implies  that  if  ABAR[I,  NK(I) ] < 0,  then 
XN0DEl1(T1)  - XNODEU(0I)  - TIME  (1 ,1)  >.  DEL  > LI  - L. 

Hence , 

XN0DEl1(Tt)  - XNODEl1(Ot)  - (LI  - L)  > TIME(I,  1) 
and 

XNODEl(Ti)  - XN0DEL(0I)  > TIME (I , 1). 

If  ABAR[I,  NK(I)]  * 0,  then  T^  would  have  been  labeled  from  0^.  Since 
feasibility  is  satisfied  at  LAMBDA  =*  LI,  it  follows  that 

TIME(I,  1)  + XNODEl1(0];)  - XN0DEL1(Tx)  < 0, 

and  consequently,  since 

ABARfl,  NK(I)  ] = TIME(I,1)  + XNODF^Oj)  - XNODE^CIj) , 

ABAR[I,  NK(I)]  cannot  be  positive. 

(lii)  Suppose  Oj  and  T^  are  both  unlabeled.  Then 

XN0DF.l(Ti)  - XNODEl(0i)  - XNODEL1(TI)  - (LI  - L)  - (XNODE^Oj)  - (LI  - L)  ] 

- XNODEL1(T1)  - XNODEL1(01) 

> TIME (I,  1) 

since  the  solution  at  LI  la  feasible. 


> 

. 
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(iv)  Suppose  0_  is  unlabeled  and  T is  labeled.  Then 


XNODEl(Ti)  - XNODE^Oj)  - XNODE^ttj)  - [XNODE^Oj)  - (LI  - L)  ] 


XNODE. 


XNODE. 


Since  feasibility  is  satisfied  at  LAMBDA  » LI 


XNODEl1(Ti)  - XNODEL1(Ox)  >_  TIME (1,1) 


and  trivially 


XNODE  (T)  - XNODE. -(0_)  + (LI  - L)  >_TIME(I,1) 


XNODE, (T)  - XNODE.  (0T)  >.TIME(I,  1) 


Now,  we  have  just  shown  that  the  XNODE  (K)'s  and  XACT  (I)'s  are 


feasible.  Lemma  6 implies  that  the  F(I,J)'s  are  always  kept  feasible. 
Lemma  4 Implies  that  the  optimality  properties  (2.42)  and  (2.43)  are 
satisfied  fbfLltfti4iferXil6DE^!(kjl',s  dfid'’¥(l,J) 's;  so  that  by  Lemma  7 these 
XNODE^(K)'i  'Sisrriatiaf^  fcomplementar$r L slackness.  Since 

w^  'flaVS  sti  6wii ti ri  £ eafiib  1 i i ii  $ Ari^  ^m^lementary  slackness  are  satis- 
fied, Theorem  3 imiikies^fifla^tfi^XNdDL9^)  V’and  xlfe'  (I) ^s^for 

Li  Li 

LI  > L >.  L2  are  optimal.  QED.  lBrno* 

, gbon  A nia  sdJ  isdmun  orIT  :fS~0$  .io') 

4.  A COMPUTER  IMPLEMENTATION  . 


A-,^qmpvjk«t  pir^gramt  implementing i!th*  improve d prdjeat  8 cfibduling<  Algorithm 


described  in  Chapter  2 is  available.  TheAtUMftLd  imput  to  the  program  is 
(a)  ^n-acyclfc  project  network  with ^OneAoUroe  and  one  sink,  and 
(^uq^ijqplJ^ti^^ftjaeeijvifetiClompifcettiMWAtliines  and  their  associated  costs 


The  program' 8 output  for  each  feasible  project  deadline  time  consists  mainly  of 

(a)  the  optimal  activity  completion  times  and  costs,  and 

(b)  the  total  project  cost. 

Optional  output  may  Include  node  labels,  optimal  node  times  for  each  project 
deadline  time,  and  dual  variables  (flows). 

Incorporated  in  this  program  is  the  option  to  have  the  minimum  project  cost 
and  corresponding  optimal  activity  completion  times  determined  for  only  one 
specific  project  deadline  time. 

A listing  of  the  computer  program  is  given  in  the  appendix.  The  flowchart 
for  this  program  is  given  in  Chapter  2,  Section  4,  pages  27-29. 


4.1.  Specific  Input  Instructions 

Card  1.  Col.  1-4:  The  number  of  nodes  in  the  network. 

Format  (14). 

Col.  6-9:  The  number  of  activities  in  the  network. 

Format  (14) . 

Col.  11:  TEST1  ■ 0 print  the  input  data, 

■ 1 do  not  print  the  input  data. 

Col.  13:  TEST2  * 9 print  the  Intermediate  output, 

■ 1 do  not  print  the  intermediate  output. 

Col.  15-18:  The  number  of  the  source  node, 

Format  (14). 

Col.  20-23;  The  number  of  the  sink  node. 

Format  (14). 

Col.  25:  TEST3  ■ 0 do  not  wish  to  specify  a single  value  for 

LAMBDA, 

■ 1 do  wish  to  specify  a single  value  for 


LAMBDA  and  print  the  intermediate  output 
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■ 2 do  wish  to  specify  a single  value  for 

LAMBDA  but  do  not  print  the  Intermediate  output. 
For  each  activity  I one  set  of  3 - 5 cards: 

Card  1.  Col.  1 - 4:  0^  * the  number  of  the  origin  node. 

Format  (14) . 

Col.  6 - 9:  Tj  ■ the  number  of  the  terminal  node. 

Format  (14). 

Col.  11-12:  NK(I)  * the  number  of  activity  completion  times  and  costs 

that  are  read  in  (jc  11)  , 

Format  (12). 

Card(s)  2-3.  Format  (8110):  TIME(I,1) TIME[I,  NK(I) ] = the  activity 

completion  times  in  increasing  order  (8 
on  Card  2 , 3 on  Card  3 if  needed) . 

Card(s)  4-5.  Format  (8110):  C0ST(I,1) , . . . ,C0ST[I ,NK(I) ] = the  cost  associated 

with  each  activity  completion  time  (8  on  Card  4, 

3 on  Card  5 if  needed) . 

The  next  card  is  present  only  if  TEST  3 = 1 or  2. 

Last  Card.  Col.  1 -10 Specified  project  deadline  time, 

LAMBDA,  Format  (110). 

The  nodes  and  activities  may  be  numbered  in  any  order.  The  current  dimensions 
will  allow  3000  nodes,  3000  activities,  and  at  most  11  different  completion  times 


and  costs. 


*T: 


i 
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4.2.  An  Example 

The  program's  input  and  output  are  illustrated  in  terms  of  the  example  net- 
work in  Figure  10.  The  input  data  are  found  in  Table  1.  As  an  example,  the  acti- 
vity cost  curve  for  activity  7 is  illustrated  in  Figure  11. 

A listing  of  the  computer  input  is  given  in  Figure  12.  The  optimal  pro- 
ject cost  curve  determined  by  the  algorithm  is  plotted  in  Figure  13.  The  opti- 
mal activity  durations  for  two  values  of  the  project  deadline  time,  LAMBDA, 
are  given  in  Table  2.  The  actual  computer  output  is  given  in  Figure  14. 


* 

1 


time 

cost 


Table  2:  OPTIMAL  PROJECT  SCHEDULES 
FOR  TWO  SPECIFIED  DEADLINE  TIMES 


Activity  # 

(I) 


Project  Deadline  Time 
LAMBDA  = 40 

Project  Deadline  Time 
LAMBDA  = 44 

Activity 

Duration 

Time 

XACT(I) 

Activity 

Cost 

Activity 

Duration 

Time 

XACT(I) 

Activity 

Cost 

68 


Figure  14: 


Computer  Output  for 


Example 


THE 

NUMBER 

OF  NODES  IS 

6. 

THE 

NUMBER 

OF  ACTIVITIES 

IS 

THE 

SOURCE 

NODE  IS  NUMBERED 

** 

NODES 

: 44 

10. 

1 ANO  THE  SINK  NODE  IS  NUMOEMEO 


6 . 


K 


2 3 4 


S 


6 


INITIAL  XNCCEIK)  0 


A IS 


16  JO 


49 


| 


| | 


*•  activities:  •• 


1 

XACT 

OHIO 

TERM 

J 

TIME 

COST 

C 

ABAR 

t 

4 

t 

2 

1 

2 

a 

0.20000E 

01 

0 

2 

4 

4 

-2 

2 

IS 

l 

3 

1 

7 

23 

0 • 30000E 

01 

0 

2 

12 

8 

0.20000E 

01 

-3 

3 

IS 

2 

>8 

3 

16 

1 

4 

1 

4 

27 

0.30000E 

01 

0 

2 

8 

IS 

0.20000E 

01 

-4 

3 

12 

7 

0. I0000E 

01 

-8 

4 

16 

3 

-12 

4 

0 

3 

4 

1 

0 

0 

O 

• 

o 

-1 

2 

0 

0 

-1 

5 

22 

2 

S 

1 

20 

a 

0.40000E 

01 

-4 

2 

21 

4 

0.3000CE 

01 

-S 

3 

22 

1 

-6 

6 

IS 

3 

S 

1 

S 

28 

0.30C00E 

01 

0 

2 

10 

13 

0.20000E 

01 

-s 

3 

IS 

3 

-10 

7 

26 

3 

6 

1 

23 

10 

0.40000E 

01 

-8 

2 

24 

6 

0.20000E 

01 

-9 

3 

2S 

4 

0. 10000E 

01 

-10 

4 

26 

3 

-11 

8 

2 S 

4 

6 

1 

23 

8 

0.20000E 

01 

-8 

2 

2S 

4 

-10 

9 

19 

S 

6 

1 

16 

12 

O.SOOOOE 

01 

0 

2 

17 

7 

0.20000E 

01 

-2 

3 

19 

3 

-3 

10 

6 

2 

3 

1 

6 

4 

O 

• 

o 

-5 

2 

6 

4 

-5 

ENTIRE 

PROJECT  COST 

CURVE 

IS  GOI 

NG  TO 

OE  DETERMINED. 

LAMHDA  a PROJECT  CCMFLETICN  TINE 
1 Hfc  STARTING  VALLE  OF  LAMBDA  IS  49. 

THE  CORRESPONDING  TOTAL  PROJECT  COST  is  0.27000E  cz. 

TmE  SOURCE  HAS  a VALUE  OF  ZERO  ANO  IS  ASSIGNEO  THE  LABEL  <-.-.-»tNF). 


ITERATION 

* 

number 

1 ••• 

the 

SINK  has  NOT  1 

SEEN  REACHEO  WITH  INFINITE  < 

CAPACITY  - CONTINUE 

WITH 

THE 

LABELING 

PROCESS. 

THE 

NOOES  THAT  HAVE  SEEN  LABELED  MILL  RETAIN 

THAT  LAUEL  FOR  THE 

REMAINDER 

OF  THE 

ITERATION. 

THE 

NODE  2 HAS 

THE  LABEL 

( 1.  1.  0. 

0.20000E  Oil. 

THE 

NOOE  1 HAS 

THE  LABEL 

( l.  1.  0. 

0.20000E  01). 

the 

NODE  4 HAS 

THE  LABEL 

< l.  1.  0. 

O.IOOOOE  0|). 

THE 

NOOE  S HAS 

THE  LABEL 

I 3,  1.  0. 

0.20000E  01 >. 

THE 

NODE  6 HAS 

THE  LABEL 

( 5,  1.  0. 

0.20000E  01). 

breakthrough:  update  the  dual  variables. 

activity  «:  1 

J 

NEW  flow:  FII.J) 

1 

1 

0.0 

1 

2 

0.0 

2 

1 

0.20000E  01 

2 

2 

0.0 

2 

3 

0.0 

3 

1 

0.0 

3 

2 

0.0 

3 

3 

0.0 

3 

4 

0.0 

4 

1 

0.0 

4 

2 

0.0 

S 

1 

0.0 

5 

2 

0.0 

3 

3 

0.0 

6 

1 

0.20000E  01 

* 

6 

2 

0.0 

6 

3 

0.0 

7 

1 

0.0 

7 

2 

0.0 

7 

3 

0.0 

7 

4 

0.0 

0 

1 

0.0 

a 

2 

0.0 

9 

1 

0.20000E  01 

9 

2 

0.0 

9 

3 

0.0 

10 

1 

0.0 

10 

2 

0.0 

ITERAT ICN 

NUMBER 

2 

the 

SINK  HAS  NOT  L1ECN  REACHEO  WITH  INFINITE  CAPACITY  - CONTINUE 

WITH 

THE 

LABELING 

PROCESS. 

the 

NODES  THAT  HAVE  BEEN  LABELED  W ILL  RETAIN 

THAT  LABEL  FOR  THE 

REMAINDER 

OF  THC 

ITERATION. 

TMF 

NODE  2 HAS 

THE  LAUEL 

I 1.  t » c. 

0.20000E  01 >• 

THE 

NOOE  4 HAS 

THE  LAUEL 

< 1.  1.  0, 

O.IOOOOE  Oil. 

I 
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nonbreakthrougf:  upoatf  the  primal  variables; 
I.F.  .DETERMINE  optimal  activity  times  pop  lambda 


46. 


OELTA  {REPRESENTED  BY  "D“ ) RANGES  FROM  0 TO  3. 

LAMBDA  RANGES  FRCM  49  TO  46. 

THE  MINIMUM  CCST  PROJECT  SCHEDULE  FOR  PROJECT  OEAOLINE 
NODE  •:  K NEv  VALUE:  XNOOEIK) 


49-d: 


0 

4 

I 5-0 
16 

30-0 


6 

49-0 

PROJECT  COMPLET  ICN 

TIME  = 

49-D. 

ACTIVITY  «! 

1 

nem  value: 

XACTI 1 > 

ACTIVITY  COST 

i 

4 

0.40000E 

01 

2 

IS-0 

0 • 2000  OE 

01  ♦ 1 

3 

16 

0.30000E 

01 

4 

0 

0.0 

5 

22 

0. 1000CE 

01 

6 

IS 

0.3000CE 

01 

7 

26 

0.30000E 

01 

8 

2S 

0.40C00E 

01 

9 

19 

0.30000E 

01 

10 

6 

0.40000E 

01 

THE 

CURRENT  VALUE 

CF  THE  PROJECT  COST  IS 

0.27000E  02  4 < 

j 

NEW 

VALUES  OF 

AfiAR 

FCfi  J=l,2 

1 

i 

I j: 

l 

2 

3 

4 

S 

1 

0 

-2 

' / 

2 

3 

0 -5 

3 

0 

-4  -a 

-12 

0.20900E  OI*D) 


4 

5 

6 

7 

a 

9 

10 


-4 

-l 

0 

-8 

-S 

0 

-2 


-4 

-2 

-S 

-9 

-7 

-2 

-2 


-3 

-10 

-10 


-l  I 


1TERATICN  NUMBER 


*** 


The  SINK  HAS  NCT  BEEN  REACHED  MlTH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LABELING  PROCESS. 
THE  NODES  THAT  HAVE  BEEN  LABELED  MILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION. 


THE 

NODE 

2 

FAS 

TFE 

LAUtL 

1 

1. 

1 . 

0. 

0.2COOOE 

on 

THE 

NODE 

3 

FAS 

THE 

LABEL 

1 

I. 

2. 

0. 

0. IOOOOE 

on 

THE 

NODE 

4 

FAS 

TFE 

LABEL 

1 

1. 

1. 

0. 

0. IOOOOE 

on 

- — --- 


i iililMMlI 


| | -if  rjA.i'1 1-*  I I — — ' “ 
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NONBRFAKTHROUGf:  UPDATE  THE  PRIMAL  VARIABLES! 

I.E.*  DETERMINE  OPTIMAL  ACTIVITV  TIMES  FOR  LAMUCA 


DELTA  (REPRESENTED  BY  "D" ) RANGES  FROM  0 TO  1. 

LAMBDA  RANGES  FRCM  «6  TO  AS. 

THE  MINIMUM  CCST  PROJECT  SCHEOULE  FOR  PROJECT  DEADLINE 
NODE  •:  K NEW  VALUE:  KNOOEIK) 


PROJECT  COMPLETICN  TIME 


ACTIVITY  » : I NEW  VALUE:  XACTIII 


ACTIVITY  COST 
O.AOOOOE  01 
0.80000E  01 
C.30C0CE  01 
0.0 

O.IOOOOE  01 

0.30000E  01  ♦ I 0#  20000E  01*0) 

0.30000E  01 

O.AOOOCE  01 

0.30000E  01 

O.AOOOOE  01 


0.33000E  02  » I 0.20000E  01*0) 


THE  CURRENT  VALUE  CF  THE  PROJECT  COST  IS 


NEM  VALUES  OF  Ae AR  FOR  J= I • 2 


ITERATION  NUMBER 


T MC  SINK  HAS  NCT  been  REACHEC  WITH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LAOELING  PROCESS. 
THE  NODES  THAT  HAVE  BEEN  LABELED  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  CF  THE  ITERATION 


0.20000E  01) 


2 FAS  ThE  LABEL  I I 


THE  NOOE 


O.IOOOOE  01) 


THE  NODE 


O.IOOCOE  01) 


A FAS  TFE  LABEL  I I 


THE  NODE 


0.20000E  01) 


5 FAS  THE  LABEL  4 2 


THE  NOOE 


NONUREAKTHPOUGF  : upcate  the  primal  variables; 

i.c.  ; determine  optimal  activity  times  FOR  LAMOOA  - 43. 


DELTA  < REPRESENTED  BV  “D" I RANGES  FROM  0 TO  2. 

LAMEOA  RANGES  FRCM  45  TC  43. 

THF  MINIMUM  CCST  PROJECT  SCHEOULE  FOR  PROJECT  DEADLINE  m *5-0 : 

NODE  »;  K NEW  value:  anode  IK ) 

1 0 

2 4 

3 12 

4 16 

5 26 

6 45-0 


THF  SINN  HAS  ACT  BEEN  REACHEO  KITH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LABELING  PROCESS. 
THF  NOOES  THAT  HAVE  BEEN  LABELED  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION. 


THE 

NOOE 

2 

FAS 

TFE 

LABEL 

I 

1. 

1. 

0. 

0 • 20  0 COG 

on 

THC 

NOOE 

3 

FAS 

THE 

LABEL 

1 

1. 

2. 

0. 

0. IOOOOG 

01) 

the 

NOOE 

4 

FAS 

THE 

LABEL 

1 

1. 

1. 

0. 

0. IGOCOE 

on 

THE 

NODE 

S 

FAS 

THE 

LABEL 

( 

2. 

I. 

0. 

o.2ocoor 

01  ) 

THE  NODE  6 HAS  THE  LAUEL  ( 


breakthrough:  update  the  dual  variables. 


0.20000E  01), 


activity  ■:  1 

i 

1 

2 
2 
2 
3 
3 
3 
3 
A 
A 
S 
S 

5 

6 
6 
6 
7 
7 
7 
7 

a 

a 

9 

9 

9 

10 

10 


new  flow:  fi i. j) 
0.2C000E  01 
0.0 

0.20000E  0| 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

0.20000E  01 
0.0 
0.0 

0 . 200 COE  01 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

0.20000E  01 
0.20000E  01 
0.0 
0.0 
0.0 


iteration  number 


the  sink  has  act  been  reached  with  infinite  capacity  - continue  with  the  labeling  process. 

THE  NODES  THAI  HAVE  BEEN  LABELEO  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION. 


THE  NODE 


THF  node 


3 FAS  THE  LABEL  I I. 


A FAS  THE  LAUEL  I 


O.IOOOOE  01). 


O.IOOOOE  01). 


nonhre akthrolgf : update  the  primal  variables: 

I.E.  OETFRMINE  OPTIMAL  ACTIVITY  TIMES  FOR  LAMBDA  = 


PELT  A I RFPRES  ENTEC  BY  "0")  RANGES  FROM  0 TC  2. 

lambda  ranges  frcm  aj  to  ai. 

THE  minimum  CCST  PROJECT  SCHEDULE  FOR  PROJECT  OEAOLINE 
NODE  •:  K NEW  value:  XNCOCIK) 

1 0 

2 A-O 

3 12 

A 16 

3 26-0 


74 


V > 


t-LLi 


43-D 


PROJECT  COPPLETICN  T IPE  = 


0-0. 


ACTIVITY  »:  I 

NEP  VALUE)  XACTII) 

ACTIVITY  COST 

1 • 

4-0 

0.40000E  01 

2 

12 

0.80C00E  01 

3 

16 

0.300C0E  01 

4 

0 

0.0 

5 

22 

0.1000CE  01 

6 

14-0 

0.50000E  01 

7 

26 

0.30000E  01 

a 

25 

0.40C0CE  31 

9 

17 

0 . 7000  CE  01 

10 

6 

0.40030E  01 

♦ ( 0.20000E  01*0) 


♦ { 0 . 20000E  01*01 


THE  CURRENT  VALUE  CF  THE  PHCJECT  COST  IS 


0.39000E  02  * ( 0.40000E  01*0). 


NEP  VALUES  OF  A8AR  FCR  J=l,2. 


1 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


j: 


«nk< i ) 
3 


2 

3 

0 

-4 

0 


0 

0 

-4 

-4 

-l 


-5 

-8 


-12 


-3 

0 

2 

-4 


-4 

-2 

0 

-4 


-2 

-7 

-5 


-1 


**» 


ITERATION  number  7 *** 

THE  NODE  2 FAS  TFE  LABEL  I 1.  2.0. INF). 


THF  SINK  HAS  NCT  BEEN  REACHEO  PITH  INFINITE  CAPACITY  - CONTINUE  PITH  THE  LABELING  PROCESS. 
THE  NODES  THAT  HAVE  BEEN  LABELEO  PILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION, 


THE  NODE 

3 

FAS 

THE 

LABEL 

( 

1.  2.  0. 

0.  IOOCOE 

01  ) 

THE  NODE 

4 

FAS 

THE 

LABEL 

1 

1.  I.  0. 

0.10000E 

01  > 

THE  NOOE 

5 

FAS 

THE 

LABEL 

1 

2.  1,  0, 

0. 10000E 

01) 

THE.  NOOE 

€ 

FAS 

TFE 

LAOEL 

( 

4,  1,  0. 

0.10000E 

01  > 

breakthrough: 

1 . ri 

UPOATE 

THE  OUAL 

variables. 

ACTIVITY 

a: 

, 

J 

Nfc  W 

flop:  fii.ji 

1 

1 

0.20000E  01 

1 

2 

0.0 

2 

1 

0.20000E  01 

2 

2 

0.0 

2 

J 

0.0 

75 


O.IOOOOE  01 
0.0 
0.0 
0.0 
0.0 
0.0 

0 • 20000E  01 
0.0 
0.0 

0.20000E  01 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

O.IOOCOE  Ot 

0.0 

0.20000E  01 
0.20000E  01 
0.0 
0.0 
0.0 


***  ITERATION  NUMBER  8 *** 

THE  NODE  2 HAS  THE  LABEL  ( t.  2.0. INF). 


THE  SINN  HAS  NCT  OEEN  REACHED  WITH  INFINITE  CAPACITY  - CONTINUE  » l TH  THE  LAbfcdNG  PROCESS. 
THE  NODES  THAT  HAVE  BEEN  LABELED  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION 


THE  NODE 


THE  NODE 


THE  NODE 


3 PAS  THE  LABEL  I 


S HAS  THE  LABEL  I 2. 


6 HAS  THE  LABEL  I 


breakthrough:  update  the  dual  variables. 


ACTIVITY  »:  I 


NEW  FI ow:  FI  I , J) 
O.20000E  01 
0.0 

0.20000E  01 
0.0 
0.0 

O.IOOOOE  01 
0.0 
0.0 
0.0 
0.0 
0.0 

0. JOCOCE  01 
0.0 
0.0 

0 . 200 COE  01 


O.IOOOOE  01). 


o.ioocoe  on. 


o.iooooe  on. 


'6 


9 

a 

9 

9 

9 

10 

10 


0.  10000E  01 

0.0 

0.20000E  01 

O.JOOCCE  01 

0.0 

0.0 

0.0 


ITERATICN  NUMBER 


The  nooe 


2 NAS  The  LABEL  ( 


• •• 

l.  2,0. INF). 


the  sink  has  nct  been  reached  with  infinite  capacity  - continue  with  the  labeling  process, 
the  nodes  that  have  been  labeleo  will  retain  that  label  for  the  remainder  of  the  iteration. 


THE  NODE 


3 HAS  THE  LABEL  I 


0.10000E  01). 


nonbreakthrolgh:  upcate  the  primal  variables; 
i.e.  determine  optimal  activity  times  for  lamboa  = 


40. 


i 


OFl  TA  (REPRESENTED  BY  "D" ) RANGES 
LAMBDA  RANGES  EBCM  A I TO 

the  minimum  cost  project  schedule  for 
nooe  «:  k new  value:  xngdeik) 


FRCM  0 TO 
40. 

project 


deadline 


PROJECT  COMPLET ICN  TIME 


0 

2 

12 

16-0 

24-D 

41-0 


41-0. 


4 i-d: 


ACTIVITY  »:  I 

NEW  value:  XACTI 1 ) 

ACTIVITY  COST 

1 

2 

o.boooce 

01 

2 

1 2 

0.8000CE 

01 

3 

16-0 

0.30000E 

01 

♦ 

( 

0. 10000E 

01*0) 

4 

0 

0.0 

5 

22-0 

0.  I0C00E 

01 

♦ 

( 

0. 30000E 

01*0) 

6 

12-0 

0 . 9000  OE 

01 

♦ 

I 

0.2000CE 

0 1 *D  ) 

7 

26 

C . 30000E 

01 

a 

25 

0.40COCE 

01 

9 

17 

0.70C0CE 

01 

10 

6 

0 • 4000  OE 

0 1 

THf  CURRENT  VALUE  OF  THE  PROJECT  COST  IS 


0.47000F  02  4 I 0.60000E  01*D>. 


HEW  VALUES  OF  AEAR  FCR  J= I . 2 ■ 


j: 


,NK( I ) 
3 


* .kitr  A A 


THE  SINK  HAS  NCT  BEEN  REACHED  KITH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LABELING  PROCESS. 
THE  NUDES  THAT  HAVE  BEEN  LAUELEO  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION. 


THE  NODE  3 FAS  THE  LABEL  I 1, 

THE  NODE  5 FAS  THE  LABEL  ( 2. 


0,  0 • 1 0 OCOE  01). 

0.  0. IOOOOE  01). 


nonbreakthrolgh:  upcate  the  primal  variables; 

l .E.  DETERMINE  OPTIMAL  ACTIVITY  TIMES  FUR  LAMBOA  = 


DELTA  I REPRESENTED  BY  "l>»  ) RANGES  FRCM  0 TC  l. 

I AMHOA  RANGES  FRCM  AO  TO  39. 

THE  MINIMUM  CCST  PROJECT  SCHECULE  FOR  PROJECT  DEAOLINE  = 

nooe  »;  k new  value:  anode (k) 

1 0 

2 2 

3 12 

A I 5-0 

5 23 

6 AO-D 


39. 


ao-o: 


PROJECT  COMPLEIICN  TIME 


AO-O. 


ACT  ! V CTY 

• : i 

new  value: 

XACT< I ) 

ACTIVITY  COST 

t 

2 

0.8000CE 

0 1 

2 

1 2 

0 • 00  OC  OE 

01 

3 

is- a 

C .40000E 

01  ♦ I 

0. IOOOOE  01*0) 

4 

0 

0.0 

5 

2 1 

0*4000  OE 

01 

t 

1 1 

0.11C0CE 

02 

i 

26 

C.  iOOOOE 

01 

8 

2b 

O.AOOOOE 

0 1 

9 

1 7-0 

0. 70C  )CE 

01  ♦ t 

O.SOOOOE  Ol*D> 

10 

6 

O.ACOOOE 

0 I 

THE 

Current 

VALUE 

CF  THE  PROJECT  COST  IS 

0.S3000F  02  ♦ ( 

0.60000E  01*0). 

Nf.  « 

VALUES 

OF  AU*R 

FOR  J=|,2, 

i j: 

i 

2 

4 

s 

6 7 

i 

2 

0 

***  iteration  number  u *** 

THE  NODE  2 FAS  THE  LAI3EL  I 1,  2.0.INFI. 

THE  SINK  HAS  NCT  3EEN  REACHEO  KITH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LABELING  PROCESS. 
The  NODES  THAT  HAVE  OEEN  LABELED  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINOER  CF  THE  ITERATION. 


THE  NODE  3 

l-AS 

THE 

LABEL 

I 1 . 2. 

0. 

0.  lOOOOE 

Oil  . 

THE  NOOE  5 

l-AS 

THE 

LABEL 

I 2 . 2 . 

0. 

0. lOOOOE 

01  I . 

• ; 

THE  NODE  6 

h AS 

THE 

LABEL 

I 5.  3. 

0. 

0. lOOOOE 

Oil. 

breakthrough: 

UPDATE 

THE  OUAL  VARIABLES. 

activity  *: 

X 

3 

NEW  flow:  FI  I 

. J) 

1 

1 

0.20000E 

01 

1 

2 

0.0 

2 

1 

0.20000E 

01 

2 

2 

0.0 

2 

3 

0.0 

3 

1 

0.  lOOOOE 

01 

3 

2 

0.0 

3 

3 

0.0 

2 

4 

0.0 

4 

1 

0.0 

4 

2 

0.0 

5 

1 

C.30000E 

0 1 

5 

2 

0. 1C000E 

0 1 

5 

3 

0.0 

• 

6 

1 

0.20C00E 

01 

6 

2 

0.0 

6 

3 

0.0 

7 

1 

0.0 

7 

2 

0.0 

? 

3 

O 

• 

o 

0.0 

0.10000E  01 

0.0 

0.20000E  01 
0.3COOOE  01 
0 . I OOOOE  01 
0.0 
0.0 


■.*+'  ^ - *"  > ' 


2 FAS  THE  LAUEL  I 1 


the  ncjoc 


THE  SINK  HAS  ACT  OEEN  REACHED  WITH  INFINITE  CAPACITY  - CONTINUE  WITH  THE  LABELING  PROCESS. 
THE  NODES  THAT  HAVE  BEEN  LABELED  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  ITERATION 


O.IOOOOE  Oil 


THE  NODE 


3 HAS  TEE  LAUEL  I 


NONQRE  AK THROUGH  I UPDATE  THE  PRIMAL  VARIABLES! 
I.E.  DETERMINE  OPTIMAL  ACTIVITY  TIMES  EUR  LAMBDA 


DELTA  IREPRFSENTEO  OY  "D"  > RANGES  FROM  0 TC  I. 

LAMHOA  RANGES  FKCM  3S  TO  3B. 

THE  MINIMUM  CCST  PROJECT  SCHEDULE  FOR  PROJECT  DEAOLINE 
NODE  »:  K NEW  VALUE!  XNODC(K) 


new  value : x act < I ) 


ACTIVITY  COST 

o.aoocoE  oi 
o.nocccc  oi 

0 • SOOOOfc  01  ♦ ( O.IOOOOE  0 1 *D ) 

0.0 

O.AOOOOE  01  ♦ I 0.40000E  01*D) 

0.11 00  OE  02  ♦ ( 0.20000E  01*01 

O.IOOOOE  01 

O.AOOOOF  01 

0.12000E  02 

O.AOCOCE  01 


THE  CURRENT  VALUE  CF  THE  PROJECT  COST  IS 


0.59000E  02  ♦ ( 0.70000E  01*0) 


NEW  VALUES  OF  AEAR  FCfi  J= 1 , 2 


ITERATION  NUMOER 


2 HAS  THE  LAHEL  I I 


ThE  NODE 


L _ 

■i  i in  ■ ' i 

. a"-.*  X--..-  j 

80 


THE  NOOE  5 HAS  TEE  LABEL  < 2.  3.0. INF). 
THE  NODE  6 HAS  THE  LABEL  ( S.  J.0.1NE). 


***** 


the  sink  was  reacheo  »ith  infinite  capacity  implying  an  infeasible  solution  to  the  primal  problem 
IF  LAMBDA  opups  eelcw  its  CURRENT  VALUE.  38 • 


I 


c ^ Appendix:  Program  Listing  COSTOOI  , 

C COSTO'02  | 

c this  program  is  designed  tc  find  the  minimum  project  cost  as  a COST002 

C FUNCTION  OF  PROJECT  DEADLINE  TIME.  CURRENT  DIMENSIONS  MILL  COSTO02  \ 

C ALLOW  A PROJECT  NETWORK  WITH  UP  TO  30CO  NODES,  3000  ACTIVITIES,  COST003 

C AND  11  LEVELS  OF  COSTS  AND  TIMES.  ALL  VARIABLES  ARE  INTEGER*?.  COST004  1 

C ( IF  ANY  VARIABLE  IS  NOT  ALREADY  IN  INTEGER  FORM,  THE  VALUES  MUST  COSTOOA 

C BE  RESCALED  - THAT  IS,  MULTIPLIED  BY  AN  APPROPR I ATF  POWER  OF  10  - CO S TOO 5 

C UNTIL  THE  VALUES  ARE  INTEGER.)  COST005 

C COST  00  6 

C** ******************************************* ************************** CO ST 00 6 
C COST 00  7 

C THE  INPUT  IS  AS  FOLLOWS  (ALL  RIGHT-JUSTIFIED)!  COSTCOT 

C CO  ST 008 

C * COLUMN  DESCRIPTION  COST008 

C CARD  1 ! 1-*  NUMBER  OF  NODES  COST009 

C 6-9  NUMBER  OF  ACTIVITIES  COST009 

C II  OPTION  TO  SUPPRESS  PRINTING  OF  INPUT  - TEST  1 COSTOIO 

C ( 0 =PR I NT * 1 =N0  PRINT)  COSTOIO 

C 13  OPTION  TO  SUPPRESS  INTERMEDIATE  OUTPUT— TEST3  COSTOl 1^ 

C I 0 = PR I NT . 1 =NO  PRINT)  COSTOll 

C 15-18  SOURCE  NCDE  C0ST012 

C 20-23  SINK  NODE  COST012 

C 25  OPTION  TC  SPECIFY  VALUE  FOR  LAMBDA  - TEST 3 C0STC13 

C 1 0=NO , 1 = YE  S AND  SEE  INTERMEDIATE  C0ST01 3' 

C OUTPUT,  2=YES  BUT  NO  INTERMEDIATE  OUTPUT ) CO  ST  0 1 A 

COST  014' 
COST 0 l 5< 
COSTOl 5 
C0ST01 6< 
COSTOl  6* 
CO  S T 0 1 7 < 
COST  017! 
CO  S T 0 1 8 ( 
COST018! 
COSTOl  9< 
CO  S T 0 1.9  i 
C05T020<: 
COST020! 
COST 021 C 
CO  ST 021 « 
COST  022C 
COST022E 
COST  023C 

C *********** *********************** ********************* ********* * »*****CO ST 02 35 

C C0ST024C 

COST024! 
COST02  50 

ABAR(I.J)  = T I ME ( I • NK ( I ) *1— J ) ♦ XNODE I OR I G ( I ) >-XNOOE( TER Ml  I ) ) COST0255 

C(I.J>  = DECREASE  IN  I TH  ACT'S  COST  PER  UNIT  FOR  J TH  TIME  COST0260 

CAP  « MINI  FLOW  REACHING  ORIGIN  NODE.  EXCESS  CAPACITY  TO  COST0265 

TERMINAL  NCOE ) COSTC27C 

COST(I.J)  * COST  OF  COMPLETING  ACTIVITY  I AT  TIME(I.J)  COST0275 

DEL  = MI N( DEL  T A 1 , DEL  T A2  I COSTC280 

OELTA1  = MI M — ABAR ( I , J)  WITH  I LABELED  AND  J UNLABELED,  COST0285 

ABAR(I.JXO)  COST  0290 

DELTA2  = MINI ABARI I , J)  WITH  I UNLABELEC  AND  J LABELEO,  COST0295 

AB  AR ( I , J ) >0 ) COST  0300 

DIPEC(J)  = DIRECTION  OF  FLOW  REACHING  NODE  J COSTOTOS 

( 0=EORW ARD , 1 =RE VERSE  ) COST03IO 


THE  FOLLOWING  CARDS  ARE  IN  SETS  OF  3' 
(ONE  SET  FOR  EACH  ACTIVITY). 


COLUMN  DESCRIPTION 

CARD  1 ! 1-4  ORIGIN  NCDE 

6-9  TERMINAL  NODE 

11-12  NUMBER  OF  ACTIVITY  COMPLETION  TIMES 

AND  COSTS  THAT  ARE  READ  IN  (<=11 ) 

CARD ( S > 2—3 ! FORMAT  8110  COMPLETION  TIMES  (8  ON  CARD  2. 

3 ON  CARC  3 IF  NEEDED ) 

CARD( S )4— 5!  FORMAT  8110  COST  ASSOCIATED  W/EACH  COMPLETION 

TIME  (8  ON  CARD  4,  3 ON  CARD  5) 


LAST  CARD  (USE  ONLY  IF  TEST3  = 1 OR  2): 

COLUMN  DESCRI FT  I CN 

1-10  SPECIFIC  VALUE  OF  LAMBDA 


DEFINITION  OF  VARIABLE 


o u u 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FLOW( It J) 
INF 


ki(i)  = 


KOUNT 
LABEL! I) 


HERE  LABELED 


OR  2 


FLOW  IN  J TH  PIECE  OF  ACTIVITY  I 
ANY  NUMBER  GREATER  THAN  MAXICAP) 

(CURRENTLY  SET  AT  (2*MAX  +1  )) 

THE  NUMBER  OF  THE  TIME-COST  PIECE  USED  IN 
LABELING  TERM* I ) FROM  GRIG* I) 

KEEPS  TRACK  OF  ORDER  IN  WHICH  NODES 

0 IF  NODE  I UNLABELEO 

1 IF  NODE  I LABELED 

SPECIFIC  VALUE  OF  LAMBDA  IF  TEST3=1 
TOTAL  NUMBER  OF  ACTIVITIES 

NUMBER  OF  DIFFERENT  TIMES  AND  COSTS  FOR  ACTIVITY  I 
TOTAL  NUMBER  OF  NODES 
ORIGIN  NODE  FOR  ACTIVITY  I 

WHERE  THE  FLOW  IS  FROM  - USED  IN  LABELING  ONLY 
PROJECT  COST  FUNCTION 
NUMBER  OF  THE  SINK  NODE 
NUMBER  OF  THE  SOURCE  NODE 
TERMINAL  NODE  FCR  ACTIVITY  I 
OPTION  TO  SUPPRESS  PRINTING  OF  INPUT 
( 0 = FR I NT # I =NO  PRINT! 

OPTION  TO  SUPPRESS  INTERMEDIATE  OUTPUT 
I 0 = FR I NT  # l = NO  PRINT) 

OPTION  TO  SPECIFY  VALUE  FOR  LAMBDA 

(0=N0t  l=YES  ANO  SEE  INTERMEDIATE  OUTPUT, 

2=YES  BUT  NO  INTERMEDIATE  OUTPUT) 

J TH  BREAKPOINT  (DURATION  TIME)  FOR  ACTIVITY  I 
ACTIVITY  DURATION  TIME 
NODE  TIME 

XNODE( ORI G ( I ) ) — XNODE*  TERM! I )),  AN  UPPER  BOUND  ON 
THE  ACTIVITY  DURATION  TIME 
INOICES 

INODE,  I TERM, I ACT, IORIG, ID  IFF, ETC. 

= NON-INDEXED  VERSIONS  OF  X NOOE ( I ) , TER M ( I ) , X ACT ( l ) , 
OR IG*  I ) , XDIFF I I ) .ETC. 


L INPUT 
NA 
NK  II) 
NN 

OR  IG  I I ) 
0RIG2I I ) 
PCOST 
SINK 
SOURCE 
TERM! I ) 
TEST! 


TEST2  = 


TEST3  = 


t i m r ( i. j) 

X ACT  I I ) 
XNODE I I ) 
XDIFFI  I ) 


I,J»K,M»N»P  = 


COST  03 1 
COST  032| 
COST  032 
COST  03  3 
COST033 
COST  034 
CO  S T 034 
COST  035 
CO  ST 035 
CO ST 036 
C0ST036 
COST  037 
C0ST037 
CO  S T 038 
COS  TO  38 
COST  039 
COST  039 
COST  0401 
COST  04  O' 
CO  S T 04  1 1 
COST  04 1 1 
COST  042 
COST  0A2f 
COST  OA  3 
CO  ST 04  3 
COSTOAAf 
C0ST044 
COST  045 
COST045J; 
COST046 
COST046I 
COST047) 
COST047 
COST 04 8f 
COST048 


C*************************** ***************************** ***************COST 04 9( 


c 

C 

C 

c 

c 

c 

c 

c 

c 

c 


COST  C491 

dimensions:  costosoi 

NN  = TOTAL  NUMBER  OF  NODES  C0ST0501 

NA  = TOTAL  NUMBER  OF  ACTIVITIES  COST051< 

MAX  = MAXI  NK(  I ) ) COST051! 

CAP! NN  ) , FLOW  I N A , M AX ) . C I NA , M AX  ) • OR  IGI NA  ) , TERM! NA ) , T I ME  I NA , MA X ) , COST 0S2I 
COST (NA  .MAX)  ,NK( MAX) , ABARI NA .MAX) ,XOIFF(NN) .XNODE (NN  ) . XACT  ( N A > . CO  ST052< 
DIRECCNN) .LABEL! NN ) ,K 1 ( NN) .ORIG2INN) . KOUNT ( NN ) . AORD I NA ) , C0ST053< 

ND( NN) ,N00 (NN) , I P ( N A ) ,CT I ME  I N A ) COSTOS3 

COST  054t 


C** ********************************************************************* COST 0543 


IMPLICIT  INTEGER*2!  A-Z  ) 

REAL*4  CAP (3000  ) .FLOW (3000 ,1 1 ) .CI3000 . i 1 ) .PCOST. INF, PCOST I • 
IKCOST , ACOST  .PNEW 

COMMON  TIME. CT I ME ♦ XNODE . OR IG. TERM , AORD .NK , NN , NA ,LM I N ,LMA X . TE ST1 
DIMENSION  CRIGI3000I  .TER  Ml 3000 ) , T I ME  I 3000 . 1 I ) .COST!  3000.1 1 > , 

INK ( 3000  ).  ABARI3000.il ) . XD IFF ( 3 000 ) .XNODE (3000)  , 

2 X ACT ( 3 000) , DIRECI 3000 ) , LABELI  3000), 

3K 1 ( 30 OC ) .ORIG2!  3000)  .KCUNTI3000)  .AORD (3000) . CT  I ME ( 3000  ) , 
4ND(3000)«  NCD ( 3000  ), IP  ( 3000) 


INPUT  DATA 


READ ( 5 , 100)  NN  ,NA , TEST  1 , TE ST2 , SOURCE .S I NK , T EST 3 


CO  ST055< 
COST055! 
COST  0S6< 
COST056 
COST057( 
COST  057) 
COST  058( 
COST  058! 
COST  059< 
COST  059! 
CO  ST 060< 
COST  06  0 
CO  ST 06 1 
CO  S T 06  1 4 


IF(T 
K3  = l 
K 2=K 
I F(K 
WRIT 
WRIT 
IFCK 
K3=K 
GO  T 
WRIT 
DO  1 
LABE 
XD  IF 
NKM1 
KK=N 
DO  9 
IF(  T 
CCI  . 
GO  T 
C<  I. 
I F(  I 

X ACT 
l F ( X 
J J- N 
A BAR 
FLOW 
CONT 
A BAR 
FLOW 
IFIT 
WRI  T 
C(  I* 
WRI  T 
CONT 
I 

DO  4 
C(  I. 
NKMI 
PCOS 
IKK* 
DO  4 
IF<K 

XI  J* 
IF  I X 
GO  T 
XI  J* 


0. 

6.220 

T1 .EC.t ) GO  TO  401 
6.150  NN.NA, SOURCE .SINK 
1=1 .NA 

.230)  ORIGt I) ,TERM(I).NK( I ) 
I ) 

.231)  (TIME ( I .J) . J=l .KN) 
,231)  I COST (I  .J)  ,J=1 « KN  ) 

UE 

ROER 

INITIAL  VALUES 
T1  .EQ. 1 ) GO  TO  193 


GT.NN)  K2=NN 
5*151)  ( K . K*K3 • K2 ) 

6,157)  (XNODE(K) ,K=K3.K2) 

GE.NN)  GO  TO  191 
1 

192 

6.152) 

1=1 ,NA 
11  = 0 

I ) =XNODE  ( OR  IG  ( I M-XNCOE  ( TERM  ( I ) ) 

<1  I 1-1 
I ) 

= 1 .NKMI 

E( I • J + l ) — T IMEC I « J ) 1 7.8.7 

=<COST( I , J) «CCST(I . J+l ) >/ (TIME( I, J+l I-TIMEC I, J ) > 

6 

*0. 

•LT.C(I.J))  I NF*C  ( I « J) 

)=XD  IFF ( I ) 

T(1 ) .LT.TIMEC I . J+l > > X ACT ( I ) =T I ME ( I , J+ 1 ) 

I >-J+l 

, J)=TIME(  I « J J ) +XC I FF ( I) 

, J > =0 
JE 

. KN I *T I ME ( I ,1  )+XDIFF(II 
« KN > = 0 

T 1 .EO.l  ) GO  TO  10 

6.153)  I , XACT ( I ) * OR I G( I ) ,TERM( I ) , ( J . T I ME ( I . J ) . COST C I , J ) , 
. ABARI  I.J  ).J=1,NKM1  I 

6,156)  KN.TIMEI I ,KN > .COST < I ,KN ) . ABAR ( I .KN ) 

JE 

AINF  + 1 • 

1 = 1 ,NA 
( I ) )=0  • 

KC I >-l 

= 0. 

K*1  .NKMI 
E.l  ) GO  TO  40 
CT1  I ) 

.GT.TIMEI I .2) ) XI J=TIME( I .2) 

41 

CT(  I ) — T I ME ( I » K ) 


IF(XlJ.LT.O)  XIJ=0 

IF!XIJ.GT.!TIME1I  ,K*1 )-TIME! I,K>  > > XI  J=TIME!  I ,K*1  )-T  IME(  I. XI 

IF! IKK .EQ. i I GO  TO  41 

IF<C( I *K)  .GT.C!I.K-1) ) GC  TC  SO 


ClJ 


PCOST I = PCO  ST l *C ! I »K  ) ♦ X I J 
GO  TO  418 
I KK= 1 

WRITE!  6,237)  1 *1 

PCOST 1= PCOST 1 «C!  I*K )*XIJ 

CONTINUE 

PCOST=PCOST*COST!  1 . 1 J+C!  I * 1 ) * T I ME  ! I , 1 J-PCCSTI 
PNEW=PCCST 


CONTINUE 
L AMBDA=LMAX 

IF  < TESTS. GE.l)  GO  TC  700 
WRITE<6.  154) 

L I NPUT  = 0 
GO  TO  9€ 

READ (5*2321  L INPUT 

IFILINPUT.LT.  LMIN)  GC  TC  705 

IFCL INPUT  .GE .LMAX)  GO  TO  704 

IF  ITEST3.E0.2)  GO  TO  724 

HR I TE (6.155)  LINPUT 

WRITE! 6*200)  LAMBDA*  PCOST 

IFC TEST2.EQ. 1 .OR.TEST3.GE.l>  GO  TO  724 

WRITE! 6,235) 

CAP! SOURCE )=  I NF 
I TER=0 

LABEL!  SOURCE  ) = 1 

IF! TEST2.E0.1  .0R.TEST3.GE. I>  GO  TO  97 

I TER  = I TER* 1 

WRI TE ! 6 • 225 ) ITER 


INITIAL  LABELING  ITERATION 


1 = 1 

J=SOURCE 
M =0 

IF  ACTIVITY  STARTS  AT  DESIGNATED  ORIGIN*  TRY  TO  LABEL, 
OTHERWISE*  CHANCE  ORIGINS. 

IF  ( OR  I G! I ) • NE • J)  GC  TO  13 
I TERM=  T ERM I I ) 

CHECK  IF  NCDE  ALREADY  LABELED  AND 
CHECK  IF  A BAR (I  »NK(  I) ) =0 • 

IF  !LABEL!  ITERM  ) .NE  .0  .OR  . AE  AR  ! I , NK  ! I M.NE.C)  GO  TO  13 
IF  NODE  NOT  ALREADY  LABELEO  ANO  A BAR ! I * NK ( I ) ) =0  * 
PROCEDE  WITH  LABELING. 

LABEL!  I TERM ) = l 
ORIG2!  I TERM) = J 
K 1 ! 1 TERM ) = NK ( I ) 

DIRECt I TERM) =0 
CAP!  ITERM  ) = INF 

IFITEST2.EC.1  .OR.TEST3  .GE.l ) GO  TO  403 
WRITE! 6,201)  I TERM.0RIG2 tl TERM) ,K1 II  TERM) 

IF  CAN  REACH  SINK,  TERMINATE  ! IMPLIES  INFEASIBLE) 

IF  ! I TERM. EQ. SINK)  GC  TO  15 
M=M+1 

KOUNT  4 M)  = ITERM 

IF  EVERY  PATH  TESTED  AND  INFINITE  FLOW  NOT  PCSSIBLE, 

GO  ON  TO  LABELING  PART!  I I). 

1=1*1 


COST092; 
COST  093 
COST093 
COST  094' 
CO ST 094 
C0ST095I 
COST  095' 
COST  096< 
COST096! 
CO  ST 09 7 1 
CO  ST  097! 
COST  098' 
CO ST 09 8! 


COST  C99! 


C0ST099'| 
COST  10o! 
COST  100^ 
COST  10 1 * 
COST  10  1 * 
COST  102! 
COST  102!! 
CO  ST  1 0 3! 
COST  103^ 
COST  1 04  li 
COST  10  45 
COST 105d 
COST  1 05? 
COST  106! 
COST106S 
COST  107* 
COST  10  7‘ 
COST 108! 
COST  108! 


COST  10  9! 


CO  S T 1 09! 
CO  S T 1 1 0 ( 
COSTllo! 
COST  1 1 1 1 j 
COST  1 1 1 ^ 1 
COST  112*  j 
COST  11  2! 

COST 1 1 3(1 
COST  11  35  j 
COST  1140 
COST  1143  \ 

cost  its!  ; 

COST  115^ 
COST  1 16^ 
COST  1 16! 
COST  117! 
COST  1175 
COST  1 1 8<i 
CO  ST  1 1 8! 
COST  11  9( 
COST  1 19* 
COST  1201 
COST  120! 
COST  1 2 1 1 
COST  12  1! 
COST  12  2! 
COST  122  j 


rt  o n no 


: 1 


• -86- 

IF  ( I ,GT  .NA  ) GO  TO  II 
GO  TO  14 

C CHANGE  DESIGNATED  ORIGINS. 

II  IF  C J «E0. SOURCE)  P= 1 
IF(P.GT.M)  GO  TC  16 

IF  ALL  LABELED  NODES  HAVE  BEEN  SCANNED  AND  NC  NEW  NODES 
HAVE  BEEN  LABELED.  GO  CN  TO  LABELING  PART  (II). 

J— KOUN  TIP) 

P=P  + l 
1=1 

GO  TO  14 

15  I F< TEST3 • GE  • 1 ) GO  TO  404 
WRI TE ( 6 .202)  LAMBDA 

404  GO  TO  999 

16  IF(TEST2,EC.l  .OR. TEST 3 .C E * 1 ) GO  TO  405 
WRI TE( 6,203) 

NEXT  LABELING  PROCEDURE 

405  1=1 
J=  SOURCE 

C AGAIN.  CHECK  ALL  CONDITIONS  FOR  LABELING 

C IE.  CHECK  IF  NODE  IS  ALREAOY  LABELED.  IF  ABAR(I.J)=0.  AND 

C IF  THE  FLOW! I. J)  IS  LESS  THAN  ITS  UPPER  BOUND. 

20  IF  !ORIG!I  I.NE.J)  GO  TO  24 
I TERM=  TERM  ( I ) 

KN=NK  ( I 1 

DO  25  K=1 , KN 

IF  (K.EO.KN)  GO  TO  27 

I F( LABEL!  I TERM ) .NE . 0 .OR . ABAR(  I.K  ) .NE  . 0 .OR  .FL 0 * ( I . K ) . GE • 

1 (C<  t ,NK(  I ) -K)-C(  I ,NK(  D-K41  ) ) )GC  TO  25 
OIREC!  I T ERM )=  0 

C CAPACITY  IS  MIN  OF  PREVIOUS  FLOW  AND  THE  EXCESS  CAPACITY 

CAP<  ITERM)=C( I ,NK( I )-K)-C( I .NK( I )-K+l  ) - FLOW(I.K) 

GO  TO  23 

27  I F ( LABEL!  I TERM) . NE . 0. OF . AB AR ( I , K ) . NE . 0 . OR  .FLOW ! I . K ) .GE  .INF) 

1 GO  TO  25 

C IF  THE  NODE  HAS  NOT  ALREACY  BEEN  LABELED.  ABAR!I.J)=0.  AND 

C THE  FLOW  IS  LESS  THAN  ITS  UPPER  BOUND.  PROCECE  WITH  THE  LABELING 

C OF  THE  NODE. 

DIRECI I TERM) =0 
CAP!  IT  ERM  ) = INF 

23  LABEL! ITERM)=1 
ORIG2!  I TERM ) = J 
K 1! ITERM )=K 

IF  ( CAP  ( I TERM  ) . GT  .CAP!  CR  IG! I ) ) ) CAP!  ITERM  ) = CAP(OR  IG(  I ).) 
IFITEST2.E0. 1 .OR. TEST3.GE  . 1 ) GO  TO  406 

WRITE (6*204)  ITERM, OR  I G2 ( ITERM  ) , K 1 ( ITERM  ).DIREC( ITERM  ) .CAP! I TERM) 
C IF  SINK  LABELED.  GO  TO  UPDATE  RRCCEOURE 

406  IF  ( ITERM. EO. SINK  I GO  TO  21 
M = M+  1 

KOUNT (M)sl TERM 

C CHECK  IF  ALL  PATHS  TRIED 

25  CONTINUE 
GO  TO  19 

24  IF! TERM!! ) .NE. J)  GO  TO  19 
t OR  I G = OR IG ( I ) 

KN=NK( I ) 

DO  26  K= 1 • KN 

IF! LAB EL!  IORIG)  .NE . 0 .OR . ABAR ! I , K ) . NE . C.OR . FLOW ( I ,K).LE.O) 

2 GO  TO  26 


COST  123 
COST  12  3 
COST  1244 
COST  124 
COST  12S 
C0STia5! 
COST  1 2 6' 
COST  1261 
COST  12  7! 
COST  127‘ 
COST  128! 
COST  128! 
COST  129< 
COS  T 1 29! 
COST130! 
COST  130! 
COST  13  1! 
COST  131! 
COST  1 32C 
COST  13  21 
COST  133! 
COST  1 3"3^ 
COST  134<] 
COST  1345 
COST  1 35? 
CO  ST  1 353 
C0ST136<; 
COST  1 363 
COST  137 
COST1 37 
COST  13  0 
COST  138 
CO  ST  1 39 
COST  1 39 
COS  T 14  0 
CO  S T 1 4 0 
COST  14  1 
CO  S T 1 4 1 
COST  14*2 
COST  142 
COST  14  3 
COST  14  3 
CO  S T 1 4 4 
COST  144 
COST  145 
COST  145 
COST 146 
COST  146 
COST  14  7 
COST 147 
COST  148 
COST  148^1 
COST  149 
COST  149 
COST  1 50 
COST  156 
COST  151 
COST  151 
COST  1520 
COST  152 
COST  153 


402 


C 

c 

c 

c 


26 

19 


18 


22 


C 

C 


407 


721 

408 


-87- 

OIREC< IOR I G )=  1 
CAP! IORIG)=FLCW( ( ,K  ) 

LABEL ( IOR I G ) = 1 
ORIG2! IOR I G ) = J 
Kill OR I G ) = K 

IF!CAP<  IORIG) .GT.CAP!  TERM!  1 ) ) ) CAP! I OR I G ) =CAP( TERM! I ) ) 

IF! TEST2.EG. 1 .GR.TEST3 .CE. 1 ) GO  TO  402 

WRITE!  6,204)  I OR  I G • OR  I G2  ( I CRI  G)  ,K1  (IORIG)  , 0 I REC  < I OR  I G ) ,C  AP  ( IOR  I G ) 
M=M  + 1 


T5 


TO  18 


KOUNT! M) =1 CRI G 
CONTINUE 
1=1*1 

IF  (I.C.T.NA)  GO 
GO  TO  20 

IF  ( J.EQ. SOURCE)  P=1 
IF  (P.GT.M)  GO  TO  22 
J=K OUNT  < P ) 

P=P+  1 
1=1 

GO  TO  20 


NONBREAK THROUGH  HAS  OCCURED.  DELTAS  ARE  FOUNO  AND  UPDATING 
MADE  IN  THE  XNCOES  AND  XACTS. 


TO  2 


DELTA1- INF41 
DELTA2*  INF  + 1 
DO  4 1=1, NA 
KN=NK( I ) 

IF  (LABEL! 0R1G (I » >.EQ. 1 .AND. LABEL! TERM! I )) .EO.O)  GO  TO  1 
A 1 IS  SET  CF  I LAEELEC  ANO  J UNLABELED, 

A2  IS  SET  OF  I UNLABELED  ANO  J LABELED. 

IF! LABEL! CRIG! I ) ) .Ed . 0 .AND .LABEL ( TERM ( I ) ) . EO . 1 ) GO 
GO  TO  4 

FINDING  0ELTA1 • S. 

DO  3 J = 1 , K N 

IF  ( ABAR! I , J)  .GE.O)  GO  TC  3 

IF  (-ABAR!  I,  J ) .LT  .DELTA1  I DELTA  1 =— AB AR  ( I , J) 

CONTINUE 
GO  TO  4 

FINDING  DELT  A2 • S 
DO  5 J= 1 , K N 

IF! ABAR! I . J) .LE .0 ) GO  TO  4 

IF  ( ABAR! I , J)  .LT.0ELTA2)  CELT  A 2=  ABAR ( I , J ) 

CONTINUE 
CONT INUE 

DEL  = MI  MDELTA1  .DELTA2) 

DEL=DELTA1 

IF  (DELTA2.LT. DEL)  0EL=DELTA2 
L AMBD A = L AMBD A — DEL 
UPDATING  THE  XNODES . 

IF! TEST2.EC.1 • OR. TEST 3 .GE • 1 J GO  TO  407 
WRITE! 6, 206 ) LAMBDA 
IF  (TEST3.F0.2)  GC  TO  721 
OEL  TA=  LAMBDA  ♦ DEL 

WRITE! 6,209)  DEL. DELT A .LAMBDA , DEL TA 

IFITEST2.E0.1. OR .TEST3.GE.1)  GO  TO  408 
WRITE! 6, 207) 

DO  80  1=1 ,NN 
I NOOE=  XNODE ( I ) 

IFlLABEL! I ) .EQ.O)  GO  TC  81 
IF! TEST2.E0.1 .0R.TEST3.GE.1)  GO  TO  409 


COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 

COST 


153 

154 

154 

155 

155 

156 

156 
157^ 

157 

158 

158 

159 

159 
160< 

160 
161 
161 
1624 
162 
163 

163 

164 

164 
1654 

165 
1664 
1664 

167 
1674 

168 
1684 
169(j 
169 
1704J 
17  01 

l 7 1 ( | 
171! 
172® 

172 

173 

1 7 341 

1741J 
1 74  31 

1750 1 
175® 
176 

176 

177 

l77l 

178® 

178 

179 

179 

180 
180 
181 
181 
182 

18 
18 
18 


IP< I ) = 2 

WRITE! 6. 21  3 I I*  ID  IFF. ACOST.C! I.FLAG2) 

GO  TO  e2 
X ACT I I 1 = I ACT 
DO  551  K = 1 . NK  Ml 
IF!K.NE.l)  GO  TO  45 
X I J=  X AC T { I ) 

IF!  XI J.GT.TIME! I ,2) I XIJ=TIME(I .2) 

GO  TO  46 

XIJ=XACT!  I )-TIME(I.K) 

IFIXIJ.LT.C)  XIJ=0 

IFIXIJ  .GT  .(TIME! I ,K  + 1 l-T IME!  I. K)  ) I XI J=TIME! I ,K  + I )-TIME<  I .K ) 
PCOST1 =PCOST 1 +C  ! I » K I *X  1 J 
CONT INUE 

KCOST  =CCST!I»1  I + CII .1  )4TIME(  I.  I I 

ACOST=KCOST-PCOST1 

PCOST=  PCOS  T + ACO  ST 

IF  CTEST3.EQ.2)  GO  TO  82 

WR I TE (6*216)  I .XACTI I ) .ACOST 

CONTINUE 

IF  ITEST3.E0.2)  GC  TC  723 

PCOSTl=PNEW 

PNEW=! FCCST-PNEWl/DEL 

WRITE! 6.224)  PCCST1  .PNEW 

PNEW=PCCST 

IF  ( TE  ST3  • NE . O.AND.LAMEDA.LE.LI NPUT I GO  TO  703 
RESET  LABELS  TO  0 AND  REFIGURE  ABARS. 

THEN  START  CVER. 

DO  87  1=1, NN 
LABEL! I >=0 
CONTI  NUE 

IFITEST2.EQ. 1. OR. TEST3.GE.il  GO  TO  420 

WRITE  16.226)  !J.J  = 1,11I 

DO  88  1=1,  NA 

NKM 1 = NK ! I )-l 

DO  500  K=1 ,KKM1 

J=NKM1  4-2-K 

ABAR! I « K ) = T IME ! I . J ) +XD IFF ! I I 
A BAR! I ,NK!  I) I = T I ME! I • 1 >+XDIFF( I ) 

IF! TEST2.EQ. I .OR.TEST3.GE. II  GO  TO  88 
NK1  =NK I I ) 

WRITE! 6,227)1 .(ABAR! I .J) ,J=1 ,NK1 ) 

CONT INUE 

IF  ! LAMBDA. LT.LMIM  GC  TC  998 
GO  TO  99 


UPDATE  THE  FLOW  AFTER  BREAKTHROUGH. 


21  IF! TEST2.EC. 1 .0R.TEST3 .GE. 1 ) GC  TO  34 
WRITE! 6.205) 

34  FLOW!  I , K 1 1 I TERM  ) )=FLOW (I ,K 1! I TERM ) )+CAP(  I TERM  I 
IF  DIREC  =0  THEN  CAP  ACDEO  TO  FLOW. 

IF  OIRIC  *1  T FEN  CAP  IS  SUBTRACTED. 

30  ITFRM*0RI62(!TFRM) 

CHECK  IF  BACK  AT  SOURCE. 

IF! ITERM. EC. SOURCE*  GO  TO  33 
FIND  WHERE  FLOW  CAME  FROM. 

I=NA 

32  1=1-1 

IF  ! OR TG!  I I.EQ.ORIG21I TERM) .AND. TERM! 1 1 .EC. ITERM)  GO  TO  31 
GO  TO  32 


COST  21 1 
COST21 j 
COST21! 
C0ST21I 

COST  2 1 1 
COST  211 
COST  21 1 
CO  ST  2 1 1 
COST2H 
COST2H 


COST  22® 
COST22(j 
COST221 
COST  221 
CO  S T 22 j 
COST22 
COST  22 
COST  22 
COST  22 
COST22 
CO  S T 22 
COST22 
C0ST22 
COST  22 
COST22 
COST22 
COST  228 
CO  S T 220 
COST  229> 
COST  229! 
COST23 
COST23 
CO  S T 231 
COST  23  I 
COST232 
COST  232 
COST  233( 
COST  233 
C0ST234 
C0ST2: 
COST  235 
COST  235 
COST  236 
COST  236 
COST  237 
COST  237 
COST  238 
COST238 
COST  239 
COST  239 
C0ST240 
CO  ST  240 
CO ST 24  L 
CO  S T 24  1 
COS  T 24  2 
COST  242 
COST  243 
C0ST24 
C0ST24* 
COST  244 
COST  24 


r»  r>  r»  n no 


P 


I i 1 


31 


33 


560 

415 

98 


998 


-90- 

CHECK  IF  DIRECTION  OF  FLO*  IS  POSITIVE  OR  NEGATIVE, 

ALSO  CHECK  IF  CAPACITY  IS  INFINITE. 

IFICAPI I TERM) .EC.INF)  GC  TC  33 
IF  I O IREC ( I TERM ) »EQ • 0 ) GC  TO  34 

FLOW ( I . Ki < ITERM) )=FLOW ( I , K 1 ( I TERM) )-CAP! I TERM) 

GO  TO  30 

RELABEL  ANO  START  OVER. 

IFITEST2  .EQ.l  .0R.TEST3 .GE . 1 ) GO  TO  415 
00  560  1=1 ,NA 
NK1=NK ( I ) 

OO  560  K=1  , NK  1 

WRITE! 6, 220 ) I ,K,FLOW! I ,K) 

OO  98  1=1 ,NN 
LABEL! I )=C 
CONTINUE 
GO  TO  99 

PROGRAM  TERMINATES  WHEN  EVENTUALLY  AN  INFINITE  FLOW  IS  ACHIEVED 
FROM  THE  SCURCE  TG  THE  SINK,  OR  WHEN  THE  VALUE  OF  LAMBDA  DROPS 
BELOW  THE  MINIMUM  LENGTH  CF  THE  NETWORK. 

IF( TES  T3.NE.0)  GO  TO  999 


LAMBDA 


705 


704 


60 

703 

706 

707 


L INPUT  ,LM  IN 

L INPUT  «LMAX 
LINPUT 


L INPUT 
LINPUT 


XACT { I ) = X ACT ( I )-D 
XACT ( I ) =XACT ( I)+0 


WRI TEI6.202) 

GO  TO  999 
WRITE! 6,233) 

GO  TO  999 
WRITE! 6,236) 

WRI TE ! 6 ,236) 

D=0 

DO  60  1=1 ,NA 
IP!  I ) = C 
GO  TO  707 
WRI TE! 6 ,234 ) 

WRITE! 6,23e> 

D=LINPUT-L AMBCA 
PCOST=0. 

DO  57  1=1, NA 

IFIIP1  I ) .EC. 1 . AND.O.GT.O ) 

IF! IP!  I ) .EO. 2 . AND.D.GT .0 ) 

PCOST 1 =0 . 

N KM  1 = N K ( I ) -1 
DO  51  K= 1 , NKM 1 
IFIK.NE.1 ) GO  TC  52 
XI J=  XAC  T ! I ) 

IFIXIJ .GT  .TIME!  1.21 ) XIJ=TIME!I 
GO  TO  53 

X I J = X ACT i I )-TIME!I,K) 

IFIXIJ.LT. 0)  XIJ=0 

IFfXIJ.GT.I TIME! I ,K+l )-TIME! I »K) )) 
PCOST l=PCOSTl+C! I ,K )*XI J 
CONTINUE 

KCOST =COST ! I , 1 )+C!I  , 1 )*TIME< 1,1) 
ACOST=KCCST-PCOST1 
WRITE! 6,21 6)  I , X ACT ! I ) ,ACOST 
PCOST=PCOST+ACOST 
WRI TE! 6 ,239 ) PCCST 
WRI TE!  6,228) 

STOP 

FORMAT!!  4 ,1X,I4,1X,I1,1X,I1.1X,I4,1X 
FORMAT! THE  NUMBER  OF  NODES  IS  • 
1TIVITIES  IS  •, 14, ••*,/, IX, ‘THE  SCURC 
2 THE  SINK  NOOE  IS  NUMBEREO  •,I4,«.*./ 
151  FORMAT ! • 0 • ,16X«'K*,7X.9!3X,I4,5X)) 


2) 


52 


53 

51 


57 

999 

100 

150 


XI J=T I ME! I ,K  + l )-T  I ME  I I ,K ) 


14,  IX, I 1 ) 

14, 

NOOE  IS  NUMBERED 

.«  -•  ,•  **  ncdes : •*• 


CO  ST245t 
C0ST246I 
COST  246! 
COST  24  7( 
COST247! 
CO  S T 24  8 1 
C0ST248! 
COST249* 
COST  249 1 
COST  250t 
CO  ST 25  0! 
C0ST251 ( 
C0ST251! 
COST  25  2! 
COS  T 25  2 ! 
COST253C 
COST2535 
CO  ST 254< 
COST  254? 
C0ST255C 
CO  ST 2551 
C0ST25‘6C  • 
CO ST 25 6? 
COST  257C 
C0ST25<7? 
C0ST25eC 
COST2585 
COST  259C 
COST259S 
C0ST260C 
COST  260? 
COST  26  1C 
COST  26 1 5 
COST262C 
C0ST2625 
COST  263C 
COST  26  3S 
COST  264  C 
CO  ST 26'4  ? 
C0ST265C 
COST  2655 
COST  2660 
COST  2665 
COST267C 
COST2675 
COST2680 
COST  2685 
COST2690 
COST  2695 
COST2700 
COST  27  05 
C0ST2710 
COST2715 
COST  27  20 
COST  27  25 
COS  T 2730 
COST  27  3 5 
OF  ACCOST27  4 0 
AND  COST  274  5 
COST2750 
COST2755 


I  52  FORMAT* .'  **  ACTIVITIES:  ** • ,// ,6X  . • I • . TX , • X ACT • , 6X , • OR  I G' * 3X . COST 27 

I • TERM • . 4 X , • J • » 6 X , 'TIME'»9X»'C0ST'» 1 4 X • *C  ' • 13X ,'ABAR'  ) COST27* 

153  FORMAT  f • • , 3X * I 4 , 3X  . 1 1 0 • 3 X • 1 4 , 3X , I 4,  ( T 39 . I 2 . 3X . 1 1 0 . 3X • I 1 0 , 3X , COST27 

IE  16.5* 3X, I 10) > COST  27 

154  FORMAT! THE  ENTIRE  PROJECT  COST  CURVE  IS  GOING  TO  BE  OETERMI NECOST27I 

ID.’)  COST  271 

155  FORMAT! THE  OPTIMAL  ACTIVITY  COMPLETION  TIMES  FOR  A SPECIFIED  COST27<; 

1 PRO JEC T OEAOLINE  TIME  = ',110,'  ARE  GOING  TO  £6  DETERMINED.*)  COST27? 

156  FORMAT!  ' * » T3S  » 12 • 3X • I 10.3X.I 1 0 . 22  X * 1 10)  C0ST28C 

157  F ORM AT  1 • 0 * • 4X.  'INITIAL  XNODE I K ) • • 3X . COST28C 

19  ! I 10.2X)  ) CO  S T 28  I 

200  FORMAT!  *0*  .‘LAMBDA  = PROJECT  COMPLETION  TIME'.//.  COST281 

COST  282 


IX. 'THE  STARTING  VALUE  OF  LAMEDA  IS 


2 .IX, 'THE  CORRESPONDING  TOTAL  PROJECT  COST  IS  • .E 1 6 . 5 . • . • ) C0ST282 

201  FORMAT! '0 '.' THE  NODE  ',14.'  HAS  THE  LABEL  ! • , I 4 , • , • . I 4 . • . 0. I NF ) . • ) COST283 

202  FORMAT! 'O'  .////. 30X,»*  * * * *',/////, IX,  COST283 

1 * THE  SINK  HAS  REACHED  WITH  INFINITE  CAPACITY  IMPLYING  ACOST284 

IN  INFEASIBLE  SOLUTION  TO  THl  PRIMAL  PROBLEM  './.20X,'IF  LAMBDA  DROCOST284 
2PS  BELOW  ITS  CURRENT  VALUE,  '.110.'.')  COST285 

203  FORMAT THE  SINK  HAS  NOT  BEEN  REACHED  WITH  INFINITE  CAPACITY  -COST285 
1 CONTINUE  WITH  THE  LABELING  PROCESS.',/  , IX. 'THE  NODES  THAT  HAVE  C0S.286 
2BEEN  LA8ELE0  WILL  RETAIN  THAT  LABEL  FOR  THE  REMAINDER  OF  THE  I TERACOST 286' 

3 T I ON . • ) COST 287 

204  FORMAT!  '0 THE  NODE  '.14,'  HAS  THE  LABEL  1 • , 14  , • , • , 14  , • , • . I 4 . • , • , COST  287 


IE16.5,'  ).'  ) 

205  FORMAT BREAKTHROUGHS  UPDATE  THE  DUAL  VARIABLES.' 
1.///.1X,'  ACTIVITY  4:  I • ,3X,« J‘ ,9X. 'NEW  FLOW:  F!I,J>'> 


COST288 

COST288' 

CQST2891 


206  FORMAT! NONBREAKTHROUGH:  UPDATE  THE  PRIMAL  V AR I ABLE S J • , / , I X , COST289 


l'I.E.  DETERMINE  OPTIMAL  ACTIVITY  TIMES  FOR  LAMBDA  = 

207  FORMAT!'  ','  NOOE  A i K'.SX.'NEW  VALUE:  XNODE ! K ) • ) 

209  FORMAT! 'DELTA  ! REPRESENTED  EY  "D"  > RANGES  FROM  0 TO1 
1 . 14, '.',/,  IX, 'LAMBDA  RANGES  FROM'.IIC,'  TC'.IIO.'.', 


COST  29  0 1 
COST  29  0! 
C0ST291 « 
COST  29 1 f 


2 /.IX,  'THE  MINIMUM  COST  PROJECT  SCHEDULE  FOR  PR0C0ST292< 

3JECT  DEADLINE  = *.I10.'-D:«>  C0ST292! 

210  FORMAT!'  • ,7X  , 1 4 , 12  X , I 10  ) COST293! 

211  FORMAT!'  • ,7X,I4,12X,I10,»-D' > COST2931 

212  FORMAT!  ‘PROJECT  COMPLETION  TIME  = • , I 1 0 , ' -D. ' . // . 1 X . C0ST294! 

1 • ACTIVITY  A:  I • ,3X , • NEW  VALUE:  X ACT ! I ) • . 9X,  • AC TI V I TY  COCOST294? 

2ST  • ) COST  29  5( 

213  FORMAT!'  « . 5X  . 1 4 , 12 X . I 1 0 , *-D ' . 9X , E 1 6 . 5 . • + I • ,E 1 3 . 5 . * *D ) • ) C0ST295? 


214  FORMAT!'  • , EX  , I 4 , 1 2 X , I 1 0 , • +D * , 9X , E 1 6 . 5 , • ♦ (',613.5.' 
216  FORMAT!'  • , 5X  , I 4 , 1 2X  , I 10 • 1 1 X ,E 1 6 . 5 ) 

220  FORMAT!'  • , 1 2 X , 1 4 , 2X , 1 2 , 7 X , E 1 6 • 5 > 

224  FORMAT!  'O'  .'  THE  CURRENT  VALUE  CF  THE  PROJECT  COST  IS 
1'  ♦ !'  ,E13.5,'*0).'  ) 


>£16.5. 


225  FORMAT! •-• ,'*** 


ITERATION  NUMBER', 16.' 


226  FORMAT! •«», 'NEW  VALUES  OF  ABAR  FCR 
1 • j:»  ,1 1 (5X.I2 ,3X>  > 

227  FORMAT!'  • . 2X  , I 4 , 7X  . 1 1 ! I 8 ,2X> ) 

228  FORMAT ( 1 HI  ) 

230  F ORM AT (I4.1X.I4.1X.I2) 

231  FORMAT (81 10) 

232  FORMAT! 110) 


J = 1 ,2,...  .NK(I)  • ,//  ,6X ,'  I'  »3X, 


COS  7295? 
COST  296< 
C0ST296? 
C0ST297C 
COST  29  7? 
COST  298  C 
COST29e? 
COST299C 
CO  S T 2995 
COST  300C 
COST  30  0 5 
COST30  1C 
COST301 5 
CO  S T 30  2 C 


233  FORMAT THE  SPECIFIED  VALUE  CF  LAMBDA ,',110.*, IS  LESS  THAN  THECOST302S 
1 MINIMUM  VALUE. '.110, '.IMPLYING  AN  INFEASIBLE  SOLUT I ON . • . // . 1 X , COST3030 


2* THE  PROBLEM  WILL  NCT  EE  WCRKEO.') 


COS  T 30  3 ? 


234  FORMAT!  • 1 '.' THE  SPECIFIED  VALUE  CF  L AMBD A , • . 1 1 0 . • , HA S BEEN  RE ACHEOCOST 30 40 

1.')  COST  3045 

235  FORMATCO'.  'THE  SOURCE  HAS  A VALUE  OF  ZERCANO  IS  ASSIGNED  THE  COST305C 


3LABEL  INF) 


COST  3055 


236  FORMAT  !•-*,' THE  SPECIFIED  VALUE  CF  LAMBDA.  ',110.',  IS  GREATER  THANCO  S T 3060 


AD-A071  498  TEXAS  A AND  H UNIV  COLLEGE  STATION  INST  OF  STATISTICS  F/G  15/5 

STATISTICAL  PERT?  AN  IMPROVED  PROJECT  SCHEDULING  ALGORITHM. (U) 

FEB  77  C S DUNN»  R L SIELKEN  N00014-76-C-0038 

UNCLASSIFIED  THEMIS-TR-55  NL 


I 

I 

I 


I 


, 


OR  EQUAL  TO  THE  MAXIMUM  VALUE  • • • 1 1 0 «•  S ' • 


2XNODE ( K ) • • S AND  XACT(I)*#S  ARE 
237  FORM  AT  C ■ • * • **  WARNING:  ACTIV 

I OST  FUNCTION; • ./.12X«* IE.  THE  C ( • • 14 • • • M) • • S ARE  NOT 
2G  • • ) 


COST 3065 
THE  ORIGINAL  COST307C 
COST  307S 

NUMBER  • .14**  HAS  A NON-CON VEX  CCOST3080 


IX. • THEREFORE . 
CPTIMAL.* ) 


238  FORMAT!  FCR  PROJECT  COMPLETION  TIME  * ».I10. • 


is:*  «//.ix. 

• ACTIVITY 


I • »3X» • NEW  value: 


IUTION 
1 

2 ST*  ) 

239  FORMAT! THE  CORRESPONDING  PROJECT  COST  IS  '.EI6.5**.'  ) 
END 

SUBROUTINE 


NON- INCREAS INCOST 308? 

COST3090 
THE  OPTIMAL  SOLCOST3095 
COST 3 100 

XACT (II* . 9X . * AC  TI V I TV 


C 

c 

c 

c 

c 

c 

c 

c 

c 


20 


10 

14 


11 


12 

13 

2 


ORDER 


THIS  SUBROUTINE  CETERMINES  THE  ORDER  IN  WHICH  TO  CONSIDER 
THE  ACTIVITIES  FOR  THE  CALCULATION  OF  THE  CRITICAL  PATH  TIME 
DIMENSIONS : 

NA=M=  THE  NUMBER  OF  ACTIVITIES  IN  THE  NETWORK 
NN=N=  THE  NUMBER  OF  NODES  IN  THE  NETWORK 

ORIGINA) * TERM IN A ) • AORD !NA ) *CT I ME I NA ) . XNODE! NN ) *ND! NN ) * NDD! NN ) . 

T (ME ( NA  .MAX)  .NK ( MAX ) 


IMPLICIT  INTEGER42! A-Z) 

COMMON  TIME • CTI ME .XNODE. OR IG • TERM. AORD .NK . NN • NA *LMI N. LMAX .TEST1 
DIMENSION  ORIG! 3000) .TERM (3000  )• AORD! 3000  ) .CT  IME! 3000  )• 

IX NODE!  3000) *ND( 3000)  .N0D!3000)  .TIMEI3000.1 1 ) .HKI3000) 

N=NN 

M=NA 

NDD! I ) = 1 
DO  5 1*2. N 

NDD! I ) * 0 <. 

DO  6 1*1 .M 


AORD! I )=0 
K*0 
MP=M+1 

DO  1 11=1, MP 

DO  20  I* 1 . N 
ND! I ) =NDD I I ) 

I 11  = 0 
IP* I 1 + 1 
DO  2 J* 1 . M 

IF (NO (ORIG! J ) ).NE.II) 

NDD! TERM! J ) )=IP 
I 11  = 1 

IF(K.EQ.O)  GO  TO  14 
DO  1 0 L=1 .K 

IF! AORD (L  ) . EQ . J ) GO  TO  11 

CONTINUE 

K*K+1 

GO  TO  13 

IF(L.EC.K)  GC  TO  2 
KM*K— 1 

DO  12  LL=L.KM 
AORD ILL) * ACRD (LL+  1 ) 

AORD! K )=J 
CONTINUE 

IF! Ill .EO.O)  GO  TC  3 


GO  TO  2 


COCOST  31 05 
COST  3 1 10 
COST  31  15 
COST  3120 
COST31  25 
COST31 30 
COST31 35 
COST  31 40 
C0ST3145 
CO  S T 3 1 5 0 
C0ST3155 
COST3160 
COST3165 
COST31 TO 
COST31 75 
CO  S T 3 1 8 0 
COST31 85 
COST3190 
C0ST3195 
COST3200 
CPST3205 
C0ST321 0 
C0ST32 1 5 
C0ST3220 
COS  T 3225 
COST 32  30 
COST  3235 
COST3240 
C0ST3245 
CO  S T 3250 
COST32$5 
COST3260 
CO  ST  3265 
COST32T0 
C0ST3275 
COST3280 
COST3285 
C0ST3290 
COST3295 
COST3300 
COST3305 
COST33 1 0 
COST3315 
COST  3320 
COST3325 
COST3330 
COST33  35 
CO  S T 334  0 


COST  3345 


L MI N=C PTI ME(  CPA  THT ) 

00  31  1=1 *KA 

NK1=NK( I ) 

CTIMEI  1 >*TIME<  I »NK 1 ) 

LMAX=CPTIME(CPATFT) 

RETURN 

END 

FUNCTION  CPTI VE(CFATHT) 


DETERMINE  THE  CRITICAL  PATH  TIME*.  CPTIME 

XNOOEdl  = EARLIEST  TIME  THAT  AN  ACTIVITY  BEGINNING  AT  NOOE  I 
CAN  COMMENCE 

DIMENSIONS  : 

NA=M=  THE  NUMBER  OF  ACTIVITIES  IN  THE  NETWORK 
NN=N=  THE  NUMBER  OF  NOCES  IN  THE  NETWORK 

ORIG(NA) • TERM ( NA 1 .AORD(NA)  .CTIME(NA) « XNODE ( NN ) • ND( NN ) « NOOINN ) • 

T IME(NA.MAX) ,NKI MAX ) 


IMPLICIT  INTEGER*  2 ( A«»Z  ) 

COMMON  TIME.CTIME.XNODE.ORIG.TERM. AORD *NK «NN • NA .LMI N .LMAX . TE ST 1 
DIMENSION  ORI  G ( 3000  ) . TER  MI3000  ) « AORD  ( 3000  I • CT  I ME  ( 3000  ) . 

1 XNODE! 3000 >,ND( 3000). NODI  3000) .TIME! 3000  * I 1)  .NK<3000) 

DO  1 1 = 1 » NK 

XNODE I I ) = 0 
DO  2 I 1 = 1 . NA 
I =AORD  I 1 1 ) 

IF! XNODE! ORIG! I ))+CTIME(  I) . GT. XNODE ( TERM!  I ) ) ) 

1 XNODE (TERM (I ) ) =XNODE ( OR IGI I ) )*CTIME(  I) 

CPTI ME = XNODE INN) 

RETURN 

END 
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